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The  absorption  of  cosmic  radio  noise  in  the  lower  Inosphere  due 
to  the  arrival  of  solar  cosmic  rays  has  been  calculated  for  several 
locations  at  the  high  magnetic  latitudes.  The  calculation  was  made 
In  three  parts;  the  volume  production  rate  of  electrons  from  ioniza¬ 
tion,  the  equilibrium  electron  density  resulting  from  the  competition 
between  production  and  recombination  processes,  and  the  absorption  of 
electro  magnetic  radiation  passing  through  the  ionized  region.  The 
calculation  was  made  for  proton  energy  spectra  in  the  form  of  power 
law  as  observed  by  balloon-borne  detectors  following  large  solar 
flares.  The  calculated  absorption  was  compared  with  observations 
obtained  by  rlometers  (relative  ionospheric  opacity  meter)  at  Point 
Barrow,  Port  Yukon,  College,  and  Farewell,  Alaska. 

For  the  solar  proton  event  of  July  14,  1959>  the  calculated  ab¬ 
sorption  based  on  balloon  data  does  not  agree  in  magnitude  nor  lati¬ 
tude  variation  with  the  observations  obtained  by  the  Alaskan  riometera, 
The  lack  of  agreement  in  absolute  magnitude  indicates  that  the  atmos¬ 
pheric  coefficients  used  in  obtaining  the  equilibrium  electron  density 
are  in  error.  In  particular,  the  colllslonal  detachment  as  well  as 
electron  attachment  coefficients  appear  to  require  considerable  re¬ 
vision.  The  lack  of  agreement  in  the  latitude  variation  substantiates 
observations  at  rocket  altitude  which  indicate  that  the  power  law 
energy  spectrum  flattens  considerably  at  proton  energies  below 
100  Mev, 
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CiJ\H'i3K  I 

Introduction 

For  years  it  has  been  known  that  a  direct  correlation  exists  between  en¬ 
hanced  solar  activity  and  changes  In  terrestrial  radio  wave  propagation.  A 
familiar  effect  Is  the  Increase  in  the  critical  frequencies  of  the  E,  FI  and  F2 
layers  and  stronger  absorption  in  the  D  layer  during  the  period  of  the  sun  spot 
maximum.  More  dramatic  changes  in  terrestrial  radio  propagation  are  often 
associated  with  flare  activity  on  the  sun.  The  Short  Wave  Fadeout  that  occurs 
on  the  stinllt  side  of  the  earth  following  a  solar  outburst  has  been  attributed 
to  an  increase  in  the  free-electron  content  of  the  absorbing  layer  of  the 
ionosphere.  This,  in  turn,  is  caused  by  the  ionizing  action  of  a  flux  of 
electromagnetic  radiation  from  the  solar  chromosphere  having  frequencies  in  the 
extreme  ultraviolet  and  soft  X-rays  bands.  Radio  transmission  "blackouts" 
occurring  in  the  polar  regions  are  often  associated  with  magnetic  storms.  Ih 
some  cases  a  transmission  blackout  may  be  caused  by  a  lowering  of  the  critlctd. 
frequency  of  the  reflecting  layers  with  the  consequent  failure  of  the  ionosphere 
to  return  signals  eaxl^ard  that  orginarlly  would  be  reflected.  loqpaired  radio 
transmission  adso  can  result  from  Increased  ionization  in  the  absorptive  reglcxis 
of  the  polar  ionosphere.  Blackouts  of  this  type  are  associated  with  magnetic 
storms  which  are  due  to  solar  flares,  the  primary  ionizing  agent  in  this  case 
being  a  swam  of  low-enexgy  particles  impinging  upon  the  upper  atmosphere.  These 
particles  are  able  to  reach  the  earth's  atmosphere  because  the  normal  Stozmer 
cutoff  mechanism  has  been  rendered  ineffective  by  the  magnetic  stom. 

Magnetic  storms  due  to  solar  flares  usually  commence  one  to  two  ds^ 
following  the  onset  time  of  the  flare;  however,  it  has  been  noticed  that  an 
Incresised  ionospheric  opacity  can  develop  well  before  the  terrestrial  magnetic 
field  shows  signs  of  any  severe  disturbance.  Because  the  geographicad 
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distribution  of  the  Ionospheric  region  of  enhanced  absorptive  power  Is  limited 
i  to  the  geomagnetle  polar  regions  In  these  events^  the  change  In  character  of  the 

upper  atmosphere  must  be  due  to  the  Ionizing  action  of  changed  particles.  These 
•  particles,  however,  must  be  much  more  energetic  than  those  which  comprise  the 

swarm  associated  with  a  magnetic  storm.  If  a  geomagnetic  cutoff  mechanism  is 
operative.  It  has  been  found  that  the  Increase  In  Ionospheric  opacity  and  con¬ 
sequent  radio  blackout  occurring  after  the  optical  slj^ting  of  the  solar  flare 
(but  well  before  the  onset  of  the  magnetic  storm)  Is  due  to  the  Ionizing  action 
of  the  Influx  of  solar  cosmic  rays  resulting  from  the  flare.  These  cosmic  rays 
are  now  known  to  consist  mainly  of  protons  and,  to  a  much  lesser  extent,  euLpha 
particles  having  energies  from  about  50  Mev  to  1  Bev.  The  role  of  solar  cosmic 
rays  In  the  production  of  radio  blackout  effects.  In  contrast  to  the  action  of 
the  other  agents  mentioned,  has  been  only  recently  recognized,  emd  extensive 
study  of  this  phenomenon  has  been  confined  to  the  pest  three  or  four  years. 

,  (1,2,3A) 

In  1953  Mltra  and  Shaln  (5)  originated  the  cosmic  radio  noise  technique 
for  obtaining  a  measure  of  the  Increase  In  opacity  of  the  upper  atmosphere 
resulting  from  various  solar  and  terrestrial  disturbances.  A  cosmic  nullo  noise 
flux  Is  Incident  continuously  upon  the  earth  from  the  galaxy  and  Is  subject  to 
continuous  monitoring.  The  cosmic  radio  noise  serves  in  the  capacity  of  a 
continuous  probe  of  the  Ionosphere  and  the  received  cosmic  noise  power  is  a 
quantitative  Indicator  of  the  electronic  content  of  its  lower  reaches.  A 
special  receiving  system  has  been  devised  to  monitor  the  cosmic  noise  level  at 
a  given  frequency.  This  Instrument,  c8Q.led  a  riometer,  will  be  discussed  more 
fully  In  a  later  section  of  this  thesis.  It  suffices  to  say  at  this  point  that 
.  the  riometer  can  provide  direct  infozmation  as  to  the  absorbing  power  of  the 


upper  atmosphere. 
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During  the  years  1957  to  ip60,  which  spanned  the  last  solar  oaxlnumj 
f  efforts  were  made  to  identify  the  energetic  particles  emitted  by  the  sun  as  a 

result  of  a  disturbance  on  Its  surface,  and  to  obtain  energy  spectra  of  the 
*  corpuscular  radiation.  This  task  Is  considerably  more  difficult  than  the  col¬ 

lection  of  cosmic  radio  noise  absorption  data  as  balloon-borne  counters  had  to 
be  aloft  during  the  time  of  the  event.  The  majority  of  particle-producing  solar 
flares  expel  relatively  low-energy  cosmic  rays.  These  stop  at  balloon  heights 
or  higher  In  the  atmosphere  and  cannot  be  detected  on  the  ground.  Moreover,  the 
Influx  Is  confined  to  the  geomagnetic  polar  regions  by  the  earth's  magnetic  field. 
Relatively  few  stations  exist  In  these  regions  whjch  are  equipped  with  the  neces¬ 
sary  facilities  to  carxy  out  this  type  of  measurement.  However,  the  BOlax  cosmic 
ray  energy  spectrum  Is  a  directly  measurable  quantity.  Ibe  Interaction  between 
the  solar  particles  and  the  earth's  atmosphere  can  be  considered  the  cause  of 
the  observed  cosmic  noise  absorption. 

Ibe  facts  set  forth  above  may  lead  one  to  pose  the  following  question: 

a 

given  the  solar  cosmic  ray  energy  spectrum,  a  quantity  directly  derivable  from 
measurement,  can  a  calculation  be  performed  that  will  yield  the  observed  absorp¬ 
tion  of  cosmic  radio  noise?  Such  a  calculation,  to  be  successful,  must  cor¬ 
rectly  predict  the  electron  production  rate  In  the  atmosphere  due  to  the  ionizing 
action  of  the  solar  cosmic  rays,  and  must  correctly  account  for  the  various 
electron  recoeiblnatlon  processes  present,  the  sections  that  follow  an  account 
Is  given  of  such  a  calculation,  the  results  obtained  are  discussed  and  compared 
with  cosmic  ray  and  cosmic  radio  noise  absorption  data  taken  during  a  solar 
proton  event. 

CHAPTER  U 

The  calculation  of  the  vertlccLL  attenuation  of  cosmic  radio  noise  at  a 
given  geographic  location  durlxig  a  solar  flare  event  was  performed  In  three 


principal  steps:  The  rate  of  production  of  free  electrons  In  the  atmosphere  due 
to  Ionizing  encounters  between  the  solar  cosmic  rays  and  the  constituent  moleeules 
of  the  atmosphere,  the  equilibrium  electron  density  resulting  from  the  Interplay 
between  the  pirlmazy  Ionizing  process  and  the  several  recoablnatlon  and  secondaury 
electron  production  mechanisms,  end  the  specific  attenuation  of  cosmic  radio 
noise  for  a  given  radio  frequency  due  to  electron-neutral  collisions  throughout 
the  atmosphere.  An  Integration  of  the  specific  attenuation  over  the  range  of 
alanospherlc  heights  gives  the  total  vertical  attenuation  of  cosmic  radio  noise. 

The  Electron  Production  Rate 

The  strong  attenuation  of  radio  waves  often  observed  following  a  solar  flare 
suggests  that  the  atmosphere  Is  rendered  conducting  In  regions  where  the  electron 
collision  frequency  Is  relatively  large.  Because  electromagnetic  radiations  from 
the  sun  capable  of  Ionizing  and  Og  are  strongly  absorbed  at  heights  of  8$  km 
or  higher  (6),  the  production  of  free  electrons  deep  In  the  atmosphere  Is  due 
almost  exclusively  to  the  action  of  penetrating  corpuscular  radiation. 

The  number  of  primary  electron-positive  Ion  pairs  formed  by  a  fast,  heavy. 
Ionizing  particle  Is  directly  proportional  to  the  energy  loss  the  particle 
suffers  traversing  a  path  in  the  atmosphere.  The  factor  of  proportlonedlty, 
according  to  Rossi  (7)^  for  protons  and  o  particles.  Is  approximately  3*6  X  10~^ 
Mev/lon  pair,  a  constant  over  an  extended  range  of  particle  energies.  !Ihe  number 
of  electron- Ion  pairs  formed  by  an  Ionizing  particle  per  cm  of  path  at  a  height 
h  Is  given  by: 

p(h)  dE 

-  (E),  (2.1) 

Q  dx 

where  d£/dx  Is  the  energy  loss  In  Mev/gm  cm~^  of  the  ionizing  particle  in  air, 
p(h)  Is  the  atmospheric  density  in  gV^n^i  and  Q  Is  the  average  energy  In  Mev 
required  to  form  an  electron-ion  pair.  Bspllclt  In  Eq.  (2.1)  Is  the  assiosption 
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that  at  all  helots  of  interest  (below  8^-90  km)  the  relative  proportion  of  the 
principal  constituents  of  the  atmosphere  remains  unchanged,  lhat  this  seems  to 
be  the  case  follows  from  rocket  measurements  made  of  the  atmosphere  below  the  D 
layer  performed  In  1952  over  VJhite  Sands,  New  Mexico,  by  groups  from  the  Uhlver> 
slty  of  Michigan  and  the  Uhiverslty  of  Huham  (u.K. ).  (8) 

For  cooqinitatlonal  purposes  the  altitude  variation  of  air  density  as  well  as 
the  pressure  and  electron  collision  frequency  were  taken  from  the  work  of  Nlcolet 
(9)*  The  values  for  p(h}  given  in  Nlcolet 's  paper  differ  very  little  from  those 
given  by  the  Rocket  Panel  (10),  some  el^ht  years  earlier. 

The  energy- loss  curve  used  In  this  thesis  was  taken  from  the  review  article 
(U)  by  Bethe  and  Ashkln.  Bailey  (12),  employing  a  different  method  for  cal¬ 
culating  the  rate  of  production  of  electrons  by  solar  cosmic  ray  particles,  used 
as  an  energy- loss  estimate  at  a  given  level  In  the  atmosphere  the  residual  par¬ 
ticle  range  divided  by  Its  energy.  Figure  1  exhibits  the  difference  between  the 
two  loss  curves  for  protons.  It  is  seen  that,  except  for  stopping  protons, 
Bailey's  energy- loss  values  represent  an  over-estimate  of  the  effect,  exceeding 
those  used  in  this  paper  by  approximately  a  factor  of  two. 

The  number  of  electron- ion  pairs  formed  in  a  unit  volume  at  a  height  h  Is 
given  by  Eq.  (2.1)  multiplied  by  a  weighting  factor,  Y(Ej9,9,h),  \rtiich  depends 
on  the  particle  spectrum  of  the  particular  event  being  analyzed,  and  Integrated 
over  the  appropriate  energy  interval  and  solid  angle: 

/  electrons  \  n  p(h)  dE 

q  (  - = -  )  -  /  / - (E)  y(E;0,9#h)  dE  dfl  .  (2.2) 

\  cm -sec  ''  Q  dx 

01 

The  quantity  Y(E;9,9,h)  appearing  in  the  integrand  represents  the  number  of 
particles  crossing  unit  area  per  second,  having  directions  of  Incidence  lying  in 
unit  solid  angle  centered  about  a  direction  fixed  by  local  zenith  and  azimuth 
angles  9,  9>  and  having  residual  energies  in  unit  energy  range  about  E  after 
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having  penetrated  to  an  atmospheric  height^  h.  If  It  can  be  assumed  that  the 
ionizing  particles,  on  the  average,  follow  approximately  rectilinear  trajectories 
in  the  atmosphere,  then  this  quantity  is  related  to  the  integral  energy  spectrum 
which,  in  turn,  can  be  derived  from  an  integral  range  spectrum  measurable  by 
balloon-borne  particle  counters. 

Such  counters,  while  ascending  into  the  upper  atmosphere  during  a  particle 
bombardment,  record  the  number  of  energetic  charged  particles  per  second  imping¬ 
ing  upon  them  at  a  given  height  h  (or  at  a  given  atmospheric  depth,  B).  From 
these  data,  one  can  ascertain  the  number  of  particles  crossing  unit  area  at  a 
height  h  per  steradlan,  per  second,  having  ranges  in  air  equal  to  or  in  excess 
of  the  corresponding  atmospheric  depth  R.  Ry  using  the  range-energy  relationship 
in  air  for  a  given  particle,  one  can  obtain  the  corresponding  integral  energy 
spectrum,  l.e.,  the  nuoiber  of  particles  crossing  unit  area  per  steradlan,  per 
second,  initially  having  energies  equal  to  or  in  excess  of  an  energy  This 
quantity  will  be  designated  j(Ej^). 

Differential  energy  spectra  derived  from  these  observations  besu:  a  simple 
relationship  to  the  welcditlng  factor,  Y(E;d,q>,h)  appearing  in  the  integrand  of 
Eq.  (2.2). 

y(E;0,<p,h)  =  N(Ej)  ,  (2.3) 

where  K(Ej^)  ■  -  dj^E^. 

E^,  the  energy  initially  possessed  by  the  particle,  is  now  considered  a 
function  of  the  particle's  residual  energy,  E,  at  hei^t  h  and  directions  6  and 
(p,  and  is  easily  obtained  for  a  given  E  and  h  from  the  range-energy  relationship 
for  the  particle  in  air  (11),  and  the  known  dependence  of  atmospheric  depth  on 
the  height.  (9) 

The  assumption  of  strict  rectllinearity  of  the  particle  trajectories  issued 
in  Eq.  (2.3)  cannot  be  Justified  completely  owing  to  the  action  of  the  earth's 
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magnetic  field  on  the  particle  motion.  However,  the  degree  to  which  Eq.  (2.3) 
represents  a  good  approximation  depends  on  the  curvature  imposed  on  the  particle 
paths  by  the  geomagnetic  field.  Accordingly,  a  separate  investigation  was  made 
of  the  actual  paths  of  fast  protons  in  the  atmosidiere  in  the  presence  of  a  mag¬ 
netic  field,  which,  in  the  small,  would  approximate  the  terrestrial  field.  The 
results  of  this  investigation,  which  Justify  the  approximation  represented  by 
(2.3},  will  be  described  in  detail  later  in  this  thesis. 

At  present  it  will  be  asstnaed  that  the  Influence  of  the  earth's  field  on 
the  particle  trajectories  can  be  neglected  and  that  the  Integral  in  Eq.  (2.2) 
can  be  rewritten: 

p(h)  dE 

-  (B)  N(E,)  dE  dfl  ,  (2.4) 

:  Q  dx 
m 

where  E^^  «  Ej^(E,h,0,<p) . 

From  balloon  measurements  during  solar  proton  events  (13,14,15),  it  has 
been  shown  that  the  integral  energy  spectra  of  solar  particles  can  be  represented 
by  a  power  law: 

J(>  Ej)  =  KEJ®  ,  (2.5) 

where  K  and  S  are  Independent  of  the  energy,  E^.  This  result  is  valid  over  the 
energy  Interval  from  80  to  4oO  Mev  for  protons.  Accordingly  N(E^)  has  a  similar 
fom: 

N(Ej)  -  ,  (2.6) 

where  S  »  7  -  1,  and  K  >  0/(7  -  1).  In  the  electron  production  rate  calculation, 
the  lower  limit  of  the  energy  integral,  E^,  is  a  function  of  the  local  polar 
angles,  0  and  9,  as  well  as  the  geomagnetic  latitude,  4.  It  is  assumed  in  this 
etviiiy  that  during  the  solar  particle  bonibardment,  at  least  until  the  onset  of 
the  magnetic  stom  that  usually  follows  such  outbursts,  the  Stunner  geomagnetic 


8 


cutoff  (as  modified  by  Quenby  and  Webber  (l6))  holds.  Thus>  on  the  centered- 
dipole  ax^roximation,  the  geomagnetic  cutoff  rigidity  is  given  by: 

M  cos^® 

p  =  I /o  p  >  (2.7 

K  [1  +  (1  -  sin  6  sin  cp 

where  P  is  the  cutoff  rigidity  for  particles  of  momentum  p  and  charge  Ze  at 
geomagnetic  latitude  6  ±b  the  local  zenith  angle,  and  q>  the  local  azimuth 
angle  increasing  clockwise  from  geomagnetic  north,  M  is  the  earth's  magnetic 
moment,  and  R  the  radial  distance  from  the  dipole  center  to  the  point  of  obser¬ 
vation.  For  d  =  0  (particles  incident  from  the  zenith  direction),  R  *  R^,  at 
the  earth's  sxirface  (6.^  X  10^  km)  and  M  =  8.06  X  10  gauss  cm^,  this  reduced 
to  the  well-known  expression: 


^  k  k  , 

P  =  ■  cos  ®  =  1^.8  cos  ®  (BV)  . 

° 

0 


(2.8) 


On  the  right-hand  side  of  Eq.  (2.7),  for  geomagnetic  latitudes  greater  than 
60  degrees,  the  minimum  value  the  square  of  the  term  in  brackets  in  the  denomln- 
ator  can  assume  is  (l  +  0.94)  =  3.8.  It  is  actually  always  very  close  to  4.0. 

To  a  good  approximation,  therefore,  at  hl^  geomagnetic  latitudes,  the  lower 
limit  of  the  energy  Integral,  £^,  is  independent  of  the  local  polar  angles  but 


depends  only  on  the  geomagnetic  latitude  so  far  as  the  Stormer  cutoff  is  con¬ 
cerned.  The  cutoff  rigidity  given  for  vertical  incidence,  Bq.  (2.8),  as  modified 
by  Quenby  and  Webber,  Is  used  throu^out  the  calculation  presented  in  this  paper. 

Another  factor  in  the  detenuination  of  E  is  the  earth's  shadow  cone.  This 

ID 

effect,  as  given  by  Kasper  (17),  la  Incorporated  into  the  solid-angle  integration 
in  the  form  of  wei^ttng  factors  for  large  zenith  angles  (d  >  70*^).  Ih  other 
words,  for  zenith  angles  greater  than  70  degrees  directions  bounded  by  "for¬ 
bidden"  angular  sectors  given  in  diagrams  in  Kasper's  paper  are  considered 
devoid  of  incident  particles.  Ihe  use  of  weighting  factors  to  account  for  the 
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earth's  shadow  cone  and  the  employment  of  the  vertical  geomagnetic  cutoff  to 
«  account  for  the  Stormer  cutoff  introduces  a  negligible  error  in  the  result  and 

avoids  the  necessity  of  an  integration  over  the  aximuth  angle.  Subject  to  the 
.  above  approximations,  the  value  of  q  for  a  solea*  proton  flvix  given  by  Eq.  (2.4) 

now  becomes: 

p(h)  pe-*'nl2  dE 

q  - -  /  a)(e)  sin  0  /  —  (E)  N(eJ  dE  d0  ,  (2.9) 

Q  ^o  'JE  (4)  dx  ^ 

m 

With  N(E^)  =  qE^  '  .  a  6uid  y  are  constants  related  to  the  incident  proton  flux, 
01(0)  is  the  shadow  cone  weighting  factor  for  protons  (a>(0)  =  2cd  for  0  S  70°  and 
<0(0)  <  2n  for  0  >  70°),  and  *  is  the  geomagnetic  latitude. 

The  integration  to  obtain  q  was  carried  out  numerically  for  three  trial 
proton  differential  energy  spectra,  with  exponents  7  =  3>  5>  and  7,  corresponding 
,  to  Integral  energy  spectra  of  the  power- law  type  with  e^qponents  smaller  by  one 

unit. 

For  each  of  these  si>ectra  q  was  computed  for  every  four  kilometers  of 
height  in  the  atmosphere,  from  32  to  100  km. 

Electron  production  rate  calculations  were  carried  out  for  geomagnetic 
latitudes  corresponding  to  geographic  locations  at  which  polar  blackout  absorp¬ 
tion  measurements  are  available.  Accordingly,  the  integration  was  perfonned  for 
the  following  Alaskan  stations:  Pt.  Barrow,  Ft.  Yukon,  College,  and  Farewell. 

The  geomagnetic  latitudes  and  cutoff  energies  for  these  locations  are  given  In 
the  following  table: 

Pt.  Barrow  ®  *  68.8°  E  =  20.8  Mev 

Ft.  Yukon  ♦  =  66.9  *  54.0  Mev 

College  ®  =  64.7°  E^  =  112.0  Mev 

Farewell  ♦  =  6l.6°  E^  =  261. 5  Mev 

For  all  three  differential  energy  spectra  the  normalizing  constant  a  was 
chosen  such  that  the  corresponding  integral  flux  of  protons  at  College,  Alaska, 
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-1  -2  -1 
was  100  particles  ster  cm  sec  . 


r  -  1 


TliuSj 

=  100, 


for  7  =  3,  S  7, 


so  that 

a  =  2.51  X  10^,  for  7=3, 

a  =  6.29  X  10^°,  for  7=5, 

a  =  1.18  X  10^^,  for  7=7- 

By  using  the  same  spectrum  for  each  of  the  four  locations  listed  above,  the  total 
vertical  attenuation  as  a  function  of  geomagnetic  latitude  was  obtained. 

The  electron  production  rate  for  7  =  3>  5,  7  and  =  112  Mev  (College, 
Alaska)  for  j(>  112  Mev)  =  100  are  shown  in  Fig.  2  as  a  function  of  atmospheric 
height.  It  can  be  seen  Immediately  that  the  production  rate  curves  are  almost 
identical  for  heights  exceeding  that  which  corresponds  to  the  stopping  depth  of 
the  protons  with  energies  near  the  Stonier  cutoff.  For  E  *  112  Mev,  the  corres- 

V 

ponding  depth  Is  10. 7  gm/cm^  which  occurs  at  a  height  of  31  kilometers. 

Recent  observations  (18)  of  heavy  nuclei,  principally  alpha  particles,  in 
solar  flare  radiation  suggest  that  the  role  of  aliha  particles  might  be  Important 
in  the  Ionospheric  absorption  events  such  as  considered  in  this  problem,  \ftille 
the  abundcmce  of  alpha  particles  relative  to  protons  Is  not  well  known  at  present, 
a  trial  calculation  was  carried  out  using  similar  rigidity  spectra  for  the  flare 
particles  aund  a  proton-alpha  ratio  of  10.  The  results  of  this  calculation  are 
shown  in  Fig.  3#  where  the  electron  production  rate  Is  given  for  a  proton  spec¬ 
trum  with  7*5  and  J^(112)  s  100,  and  an  alpha  particle  spectrum  with  7  *  ^•73> 
and  J^(119.5)  =  10.  The  two  electron  production  rates  are  shown  separately  as 
well  as  their  sum,  representing  the  total  electron  production  rate  during  a 
solar  flare  event  with  a  significant  cdpha- particle  content. 


11 


CHAPTER  III 

In  the  atoosfhere  above  the  lower  boundary  of  the  D  region  (which,  to  be 
definite,  will  be  assumed  to  lie  at  about  90  km),  the  recombination  of  electrons 
is  chiefly  radiative  and  production  and  recombination  processes  are  direct  and 
simply  competitive.  In  contrast  to  this,  the  region  below  ninety  km  is  charac¬ 
terized  by  nonradiative  recombination  processes  and  by  the  fonnation  of  negative 
molecular  ions.  This  is  due  to  the  fact  that  the  atmosphere  above  ninety  km  does 
more  than  simply  shield  the  lower  atmosphere  from  ionizing  radiation.  That  part 
of  the  solar  spectrum  that  causes  molecular  photodlssoclatlon  is  also  heavily 
absoibed  (6,8)  emd  free  electrons  created  in  the  course  of  a  corpuscular  bombard¬ 
ment  find  themselves  in  an  atmosphere  composed  almost  exclusively  of  the  molecular 
species  of  its  constituent  elements.  Ibe  ambient  temperature  of  the  atmosphere 
in  the  height  range  thirty  to  ninety  km  esdilbits  little  variation  around  a  mean 
value  of  about  250  degrees  K.  The  collision  frequency  is  quite  large  (10^®  to 
10^  sec'^  and  thermalization  of  the  nascent  free  electrons  occurs  in  a  matter  of 
microseconds.  An  environment  such  as  this,  slow  electrons  in  a  molecular  atmos¬ 
phere,  favors  nonradiative  recombination  processes. 

Mathematically,  the  temporal  behavior  of  the  electron  density  in  regions 
above  ninety  km  can  be  characterized  by  the  familiar  equation: 

dN/dt  =  q  -  ,  (3-1) 

where  q  is  the  electron  production  rate,  N  is  the  electron  density,  and  a  is 
the  recombination  rate  coefficient.  The  equilibrium  electron  density  profile 
can  be  obtained  trom  the  relatively  simple  relationship: 

N  =  (q/a)^^  .  (3.2) 

The  situation  is  not  so  simple  in  the  region  below  the  D  layer.  The  presence 
of  negative  ions  adds  a  new  dimension  to  the  overall  interaction  that  establishes 
the  equllibrliim  electron  density.  However,  a  relationship  analogous  to  Eq.  (3*2) 
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holds  in  this  region  eilso^  but  before  It  Is  presented,  those  Individual  con¬ 
tributory  processes,  believed  to  be  the  most  important,  vill  be  reviewed 
separately  along  with  the  values  assumed  for  their  associated  rate  coefficients. 
Electron  Loss  Mfechanlsms 

In  an  atmosphere  70St  Ng  and  2256  Og  with  other  trace  constituents  and  at  a 
temperature  of  2^0  degrees  K  the  following  are  considered  to  be  the  primary 
mechanisms  for  recombination: 

(a)  Dissociative  Recombination:  Rate  Coefficient  cx^.  Ihis  process  can  be 
described  by  the  equation 


(XY)  +  e  -♦X  +  Y  , 

where  X  and  Y  represent  the  constituent  atoms  of  the  positive  molecular  Ion, 

+  « 

XY  ,  and  X  represents  one  of  the  atoms  In  an  excited  state.  In  this  process  a 
free  electron  combines  with  a  positive  molecular  ion,  the  excess  energy  serving 
to  dissociate  the  molecule.  This  Is  a  relatively  fast  process,  rate  coefficients 
for  molecular  nitrogen  being  quoted  as  lying  between  10“^  to  lo"^  cmVsec. 

Nlcolet  (6)  gives  citp(Ng)  =  5  X  10“^  cxi^/bbc,  OpCOg)  =  3  X  lO”®  ct^/sec,  and 
c«jj(K0)  »  3  X  10  ^  cm V sec.  Crain  (19),  in  a  more  recent  paper,  quotes  various 
values  for  Ng  ranging  from  lo"^  to  4  x  lO’^  and  3  x  lo”^  to  lO"®  for  Og.  In 
addition  to  this,  another  process  of  unknown  rate,  namely 

»2^  +  Og  -» Ng  +  Og"^  , 

occurs  which  makes  It  difficult  to  determine  which  molecular  ion  Is  the  most 

important  in  the  dissociative  recooibinatlon  process.  Crain,  however,  states 

"7  “8 

that  a  value  of  10~ '  for  molecular  oxygen  is  more  reeusonable  than  10**  .  This, 

in  effect,  means  that  there  is  no  real  practical  difference  which  molecular  Ion 
Is  chiefly  responsible  for  this  reaction  with  electrons.  For  purposes  of  cal¬ 
culation,  Nlcolet 's  value  of  OpCNg)  »  ^  X  lO'*^  cm®/sec  has  been  assumed  for  the 


rate  of  dissociative  recombination  In  air.  Bailey,  however,  used  »  3  x  10 
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cmVsec  in  his  paper  on  cosmic  ray  induced  radio  vave  absorption.  (12) 

(Td)  Nonradiative  Attachment;  Rate  Coefficient  n*  Ihe  water  vapor  molecule 
and  oxygen  molecule  have  an  affinity  for  free  electrons.  0^  is  far  more  abundant 
than  HgO  in  the  upper  atmosphere  (I9)j  the  formation  of  0^  alone  can  be  con¬ 
sidered  the  one  iJig>ortant  process  of  this  type.  Molecular  oxygen  has  an  electron 
affinity  of  about  0.*t'  ev  (20),  but  the  attachment  process  itself  has  a  pressure- 
dependent  rate  as  has  been  shown  by  Chenin,  Rielps  Euad  Biondl  (21).  This, 
therefore,  means  that  a  third  body  is  involved  in  the  reaction  to  conserve 
momentum  and  carry  off  excess  energy.  !nie  process  can  be  described  by  the 
equation 

Og  +  e  +  (XY)  -*02"  +  (XY)  , 
and  the  rate  coefficient  quoted  by  these  authors  is 

Tj  a  1.5  X  10“^®  n(0g)  cmVsec  , 

where  ^(Og)  is  the  nuzaiber  density  of  molecular  oxygen,  this  species  being 
assumed  to  be  the  third  body  in  the  interaction.  Bailey  (12)  used  a  somewhat 
larger  rate  coefficient  for  this  process.  Although  it  is  not  explicitly  given 
in  his  paper,  he  assumed  z)  of  approximately  2  x  lo"^  0(0^)  cm^sec. 

Secondary  Electron  Production  Mechanisms 

Besides  the  direct  production  of  electrons  by  solar  cosmic  rays  \dilch 
initiates  the  secondary  processes  being  described,  there  are  two  other  Important 
processes  occurring  within  the  atmosphere  below  90  kn  which  produce  electrons. 
These  are  the  photodetachment  of  electrons  from  negative  ions  (mainly  O^**)  and 
coUlslonal  detachment  of  electrons  from  the  same  ion. 

Fhotodetacbment  of  negative  lens;  Rate  Coefficient  pS.  As  mentioned 
eeurller,  the  portion  of  the  solar  spectrum  which  penetrates  the  atmosphei-e  to 
depths  greater  than  the  D  layer  does  not  have  quanta  of  sufficient  energy  to 
photodissociate  neutral  molecules  or  form  electron- positive  ion  pairs.  However, 
the  negative  ion,  0^ ,  having  been  formed  in  the  three-body  process  mentioned 


earlier  nay  undergo  a  further  process  described  as  follows: 

Og”  +  hv  -♦  Og  +  e  . 

The  cross  section  for  this  reaction  for  various  frequencies  has  been  measured  by 
Smith,  Burch,  and  Branscosib  (22).  Ihelr  measurements  lead  to  a  photodlssoclatlon 
rate  of  Og*  In  the  upper  atmosphere  of  0.4h  ca^/sec.  This  particular  rate 
coefficient  seems  to  be  accepted  by  other  workers  (6,19),  and  In  that  sense  can 
be  said  to  be  unique. 

Colllslonal  Detachment  of  Electrons  from  0g~:  Bate  Coefficient  k.  Direct 
collisions  between  the  negative  Ion,  Og",  and  neutral  molecviles  In  the  atmosphere 
can  resxilt  In  the  reaction 

Og'  +  (XY)  -4  Og  +  (XY)  +  e  . 

The  value  of  the  rate  coefficient  of  this  process  seems  to  be  quite  uncertain. 
Branscooib  and  Bailey  (19)>  on  the  basis  of  Ionosphere  absorption  measurements, 
quote  a  value  of  k  =  2  x  lo”^^  cmVsec.  Nlcolet  and  Aiken  (6)  give  an  upper 
limit  of  a  a  4  X  lO'^^  cmVsec,  but  Hielps  and  Pack  (20),  on  the  basis  of  drift 
tube  swam  measurements,  present  a  value  of  this  rate  coefficient  three  orders 
of  magnitude  lower  than  those  quoted  above.  These  latter  authors  give  a  rate 
of  H  =  4  X  10'^  cmVsec  at  230®K. 

The  value  of  the  colllslonal  detachment  rate  given  by  Fhelps  and  Pack 
differs  qualitatively  as  well  m  quantitatively  from  the  others  quoted  In  that 
It  was  derived  from  a  direct  laboratory  measurement  rather  than  Inferred  from 
the  results  of  measurements  made  on  other  phenomena.  Yet,  uziless  there  exists 
an  important  source  of  electrons  In  the  upper  atmosphere,  operative  during  a 
proton  event,  that  has  been  completely  overlooked,  this  small  colllslonal 
detachment  rate  cannot  be  reconciled  with  the  observed  ni{(httlme  absorption  of 
cosmic  radio  noise.  The  cosmic  radio  noise  absorption  dc^- night  ratio  defined  as 

D  Toted  Daytime  Absorption  (db) 


N  Total  Nighttime  Absorption  (db) 
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is  known  from  the  College  riocieter  data  (23)  to  he  in  the  nelghhortiood  of  5. 
Accordingly,  a  calculation  was  made  wherein  the  value  of  k  was  varied  between 
the  extreme  values  quoted  above  and  d/n  as  a  function  of  flux  at  cutoff  was  cal¬ 
culated.  The  flux  values  assumed  were 

J(112  Mev)  »  100,  1000,  and  10,000  , 

1 8 

for  7=5  differential  proton  energy  spectrum.  It  was  found  that  a  =  5  x  lO' 

cmVsec  gives  d/n  *  5  for  eill  three  values  of  J.  The  value  for  the  coUlslonal 

"  l8  ^ 

detachment  rate  used  in  this  thesis  is  a  =  5  X  lO"  cm^sec. 

Mutual  Neutralization  of  Positive  and  Negative  Ions;  Rate  Coefficient 
Such  a  process  in  the  case  of  the  oxygen  molecular  ion,  can  be  described  by  the 
equation 

Og"  +  (xy)^  Og  +  (XY)  . 

This  process,  which  does  not  directly  affect  the  free-electxon  density  but  bears 
strongly  on  the  negative- ion- to-electron  ratio,  was  considered  an  important  one 
by  Bates  and  Massey  (24).  Nlcolet  (6)  quotes  a  value  of  «  lO’"^  cm^sec. 

This  value  of  is  the  one  adopted  in  this  thesis.  Bailey  (12),  however,  used 
a  value  of  of  approximately  two  orders  of  magnitude  less  than  the  above, 
namely  ■  3  X  lO"^  cm^/sec. 

In  the  following  table  the  values  of  the  various  rate  coefficients  adopted 
for  use  in  this  thesis  are  sumnmrlzed: 

Dissociative  Reconblnation  =  5  x  10*'^  cm^sec 

Ronradlatlve  Attachment  tj  «  1.5  x  lO"^®  n(0g)  cm^/sec 

Riotodetachment  from  Og'  pS  =  0.44  sec 

CoUlslonal  Detachment  from  Og'  a  >  5  x  10~^^  cm^/aec 

Ionic  Mutual  neutralization  «  10~^  cm^sec 

The  Effective  Becomibination  Coefficient 

Having  identified  the  principal  processes  by  which  free  electrons  interact 
with  the  various  particles  cOTprislng  the  atmosphere  in  the  hel^t  range  thirty 
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to  ninety  Ism,  one  can  express  the  over-all  teiqporal  variation  of  the  electron 
and  ion  denaitlea  in  terms  of  the  rates  of  production  and  loss  by  these  processes 
as  foUotfs: 

dH/dt  =  q  -  djflS*  -  iin(02)N  +  «nN"  +  pSN"  J  (3*3) 

dNVdt  =  q  -  OjjNN^  -  Oj^nV  j  (3.^) 

dN'/dt  =  qn(02)N  -  unH"  -  Oj^nV  -  pSN"  ,  (3*5) 

where  and  n"  are  the  electron,  positive  ion,  and  negative  ion  number  den¬ 

sities,  respectively;  n(02)  le  the  nuo^r  density  of  neutral  oxygen  molecules, 
n  is  the  total  neutral  particle  number  density,  and  q  is  the  production  rate  of 
electron-positive  ion  pairs  due  to  the  solar  cosmic  ray  flux.  The  various  rate 
coefficients  are  defined  in  the  table  at  the  end  of  the  previous  subsection. 

It  is  lomdlately  apparent  that  these  equations  are  consistent  with  the 
physical  requirement  of  charge  neutrality;  l.e.,  N  +  N~  -  >  0.  OSierefore, 

if  this  condition  is  e3q>llcltly  Imposed,  we  have  only  two  Independent  equations 
—the  third  being  directly  derivable  from  the  other  two. 

*4*  ^ 

!nie  quantities  N  and  N  can  be  expressed  in  terms  of  N  and  X,  the  negative- 
ion-  to-electron  ratio,  as  follows: 

n"  »  X»  ;  (3.6) 

-  H"  +  K  »  (1  +  X)N  .  (3»7) 

!n:e  Important  siiqplifying  assumiptlon  is  now  made  that  the  time  derivatives 
in  Eqs.  (3.3),  (3*^)  and  (3.^)  differ  significantly  from  zero  during  transitional 
periods  that  are  relatively  brief  compared  to  the  total  time  scale  of  the  event. 
If  the  substitutions  Implied  by  Eqs.  (3*6)  and  (3*7)  made  in  (3*3)  and  (3*^) 
and  the  assumption  of  quasi-equilibrium  is  made,  two  equations  relating  the 
unknowns,  K  and  X,  to  known  quantites  follow  directly. 

q  -  0(p(l  +  X)B^  -  Iin(0g)»  +  tenXN  +  pSXN  »  0  ; 

-ajX(i  +  X)li^  +  nn(0g)ll  -  itnXN  -  pSXU  -  0  . 


(3.8) 

(3.9) 
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Adding  these  tvo  equations,  one  obtains  an  expression  for  N,  the  electron  density, 
entirely  analogous  to  that  holding  for  the  atmospheric  region  above  ninety  km 

lEq.  (3.2)]. 

N  »  Iq/(1  +  X)(ajj  +  ,  (3.10) 

wherein  the  denominator,  (l  +  l)(o£p  +  Xa^),  can  be  identified  with  a,  the  recom¬ 
bination  coefficient.  In  Bq.  (3*2). 

The  complication  Introduced  Into  the  analysis  by  negative- ion  fozmatlon 
expresses  Itself  In  terns  of  the  unknown  value  of  X  In  Eq.  (3. 10).  Bailey,  In 
his  analysis  of  solar  cosmic  ray  induced  absorption,  obtained  an  expression  for 
X  In  texms  of  the  rate  coefficients  and  the  neutral-particle  densities  only. 

This  can  be  done  If  the  tern  a^X(l  -f  X)N^  is  negligible  conopared  to  the  other 
tenns  in  Bq.  (3*9)  •  Assuming  this  to  be  the  case,  one  inaedlately  obtains 

Xg  -  ijn(0g)/(im  +  pS)  .  (3*11) 

The  values  of  the  rate  coefficients  appearing  In  the  literature  subsequent  to  the 
publication  of  Bailey's  work  tend  to  place  this  assunqptlon  on  less  firm  ground. 

has  been  Increased  two  orders  of  magnitude  while  q  and  k  have  new  lover  values. 
That  this  approximation  may  not  be  valid  Is  particularly  true  for  nl^ttlme 
absorption  In  which  ceise  the  photodetachaaent  rate  coefficient,  pS,  vanishes. 
However,  Xg,  as  defined  by  Bq.  (3>ll);la  an  Important  parameter  In  obtaining  the 
equilibrium  negative- Ion- to-electron  ratio  by  the  method  used  In  this  thesis. 

This  present  method  requires  no  assumption  as  to  the  relative  magnitude  of  the 
terms  In  Eqs.  (3*8)  (3*9) • 

The  following  parameters  are  defined: 

P  ■  “i/®!  ' 

I  -  l/Xg  -  (an  +  pS)/ijn(02)  ; 

C  ■  otj^q/WOg))^  . 

Following  some  straightforward  but  sll^tly  lengtl^  algebraic  manipulations, 
N  can  be  eliminated  between  Eqs.  (3<3)  and  (3*9)  and  the  following  cubic  equation 
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obtained  in  X: 

X^  +  uX^  +  vX  +  w  =  0  ,  (3 ‘is) 

where 

u  -  [(^  -  2/06^  -  C]/(5^  -  C)  , 

V  -  (1  -  2pO/(i^  -  C)  , 
w  »  p/(S^  -  C)  • 

X,  which  now  appears  as  a  root  of  a  cubic  equation  (3 •12);  depends  on  q, 
the  electron  production  rate  of  the  particular  event;  as  well  as  on  the  atmos¬ 
pheric  parameters.  However;  It  Is  of  interest  to  note  the  behavior  of  the  cubic 
as  the  parameter  C  tends  to  zero. 

C  =  a^o/t^nCOg)]^  ->  0  ; 
u  ->  O  -  2X^)  , 

^ 

w  -»PX^2  ; 
and 

X^  +  0  -  2X^)X^  +  (Xj2  -  2pX^)X  +  pXj®  »  0  ; 

Identically  if  X  >  X^. 

For  various  atmospheric  heights  It  Is  found  that  the  discriminant  of  the 
cubic  Is  often  negative.  Ibis  guarantees  the  equation  has  three  real  roots;  one 
root  Is  negative;  having  no  physical  significance;  but  the  other  two  are  positive 
and  invariably  He  on  either  side  of  X^. 

X^  <  0;  but  0  <  Xg  <  Xg  <  X^  . 

Because  X^  Is  the  limiting  case  for  very  small  q;  the  root  ^Ich  represents  the 
actual  negatlve-lon-to-electron  ratio  must  be  less  than  X^;  as  q  represents  a 
source  of  electrons.  This  choice;  based  on  physical  reasoning;  Is  supported  by 
the  fact  that  for  cases  of  a  positive  discriminant  wherein  only  one  real  root 
exists;  Its  value  Is  Invariably  less  than  that  of  the  corresponding  X^. 
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Calculations  have  shown  that  for  the  values  adopted  for  the  rate  coeffici¬ 
ents  in  this  thesis,  Bailey's  approximation  for  X  holds  up  remarkably  well  for 
the  nighttime  case  as  well  as  the  daytime  for  edl  reasonable  cosmic  ray  flux 
values.  It  was  only  during  the  search  for  a  vadue  of  k  consistent  with  the 
observed  day-night  ratio  that  significant  differences  were  noted  between  Xg  and 

the  proper  root  of  the  cubic.  These  differences  were  of  an  order  of  magnitude 

-1  2  -1 

at  heights  of  thirty  to  forty  km  for  j(ll2  Mev)  i  500  particles  ster  cm  sec 
under  nighttime  conditions  (pS  =  0)  for  values  of  k  less  than  lo”^^  Ga^/aec.  For 
larger  values  of  k  (assuming,  of  course,  that  the  values  of  the  other  rate  coef¬ 
ficients  are  not  significantly  revised),  Bailey's  approximation  for  X  can  be  made 
with  confidence.  It  is  to  be  noted  further  that  the  quantitative  differences 
between  Bailey's  results  and  those  obtained  from  the  present  analysis  exist  by 
virtue  of  the  different  values  adopted  respectively  for  the  various  rate  coef¬ 
ficients.  Ihe  difference  in  method  of  obtaining  X  contributed  little  in  this 
respect,  although  the  method  vised  in  this  thesis  is  on  sounder  theoretical 
ground. 

In  the  literature  of  the  lower  ionosphere  the  quantity  (Op  +  Xa^^)  is  termed 
the  effective  reconiblnation  coefficient  rather  them  the  more  significant  quan¬ 
tity  (l  X)(aep  XOji^).  This  convention  is  adhered  to  in  this  thesis.  In  eon- 
foxmlty  with  this,  Bq.  (3*10)  is  rewritten: 

N  =  [q/(l  +  .  (3.13) 

where 

“eff  *  ^  ^1^  • 

^  Figs.  5,  6,  7»  and  6  the  quantities  X,  (l  ^)^eff'  ^  calculated  in 

this  thesis,  are  displayed  along  with  their  counteziiarts  as  ccsqputed  by  Bailey 
for  day  and  night.  Bailey's  results  are  designated  B,  whereas  those  obtained 
herein  are  designated  A.  1b  plotted  from  right  to  left  instead  of  in  the 

conventional  manner  for  the  sedse  of  clarity  of  presentation  of  the  results 
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displayed  In  Figs.  5  and  6.  The  ratio 

-  [(1  .  ' 

Where  here  refers  to  Bailey's  computed  values  for  the  negative- ion- to-electron 
ratio,  is  the  ratio  of  electron  densities  for  a  given  production  rate  6us  would 
be  obtained  using  the  rate  coefficients  of  Bailey  and  those  adopted  in  this 
thesis  respectively.  It  is  presented  for  day  emd  night  in  Fig.  9> 

The  EOLuilibrium  Electron  Density 

The  equilibrium  electron  density,  which  is  of  central  importance  to  this 
study,  is  obtained  Inmediately  from  Eq.  (3*10)  or  (3*13);  assimilng  that  the 
production  rate  and  the  negative- ion-to-electron  ratio  are  known.  Figure  4  gives 
the  daytime  and  nighttime  electron- density  profiles  due  to  a  pure  proton  flux 
[j(ll2  Mev)  -  100;  7  =  5l  and  a  proton-alpha- particle  flux  in  a  particle  ratio 
of  10  to  1  [J„(112  Mev)  =  100;  7^  =  5  and  J  (119.5  Ifev)  «  10;  7  =  4.73].  Con- 
tours  of  the  equilibrium  electron  density  as  a  function  of  geomagnetic  cutoff 
for  constant  atmospheric  heights  (daytime)  are  presented  in  Figs.  17>  I6,  and  19* 

CHAPTER  IV 

Validity  of  the  Equilibrium  Assumption 

One  of  the  most  Important  and  useful  approximations  made  in  the  development 
of  the  method  presented  in  this  thesis  for  the  determination  of  the  electron 
density  belf^t  profile  from  the  electron  production  rate  resides  In  the  assumption 
of  equilibrium.  Equilibrium  in  its  strictest  sense  Implies  simultaneity  with 
respect  to  changes  of  q  and  N.  On  the  other  hand,  these  changes  sust  take  place 
with  extreme  slowness,  as  only  In  this  manner  can  dN/dt  be  strictly  zero. 

Actually,  dH/dt  Is  not  assumed  to  actually  vanish  but  always  to  assume  values 
that  are  negligibly  small  compared  to  the  magnitudes  of  the  various  terms  repre¬ 
senting  the  several  electron  production  and  recombination  processes.  An  assump¬ 
tion  such  as  this  might  more  properly  be  termed  one  of  quasi-equilibrium.  Here, 
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as  q  changes  In  tine  at  some  finite  rate  at  a  given  height,  N  is  assumed  to  pass 
through  a  sequence  of  equlllhrlum  states  corresponding,  at  each  moment,  to  the 
particular  value  of  q  related  to  N  through  the  equation  of  eq\illlbrlum.  An 
estimate  of  the  degree  of  validity  one  can  ascribe  to  the  quasi- equilibrium 
approximation  can  be  obtained  from  the  determination  of  the  lag  time  (or  the 
amount  of  departxire  from  simultaneity)  exhibited  between  changes  in  q  and  the 
concomitant  changes  in  N*  If  the  response  time  between  N  and  q  is  small  relative 
to  some  characteristic  time  Interval  of  the  event,  changes  in  N  brought  about  by 
changes  in  q  can  be  considered  almost  simultaneous  and  quasi- equilibrium  is  a 
valid  approximation.  If,  on  the  other  hand,  the  response  time  is  relatively 
large,  a  quasl-equlllbrlum  process  cannot  be  assumed. 

The  purpose  of  this  subsection  is  to  provide  some  quantitative  infozmatlon 
concerning  the  validity  of  the  quasi-equilibrium  approximation  used  in  this  thesis. 

Two  idealized  sets  of  initial  conditions  are  proposed  that  would  result  in 
the  maxlmuBi  departures  fxtim  equilibrium  between  the  electron  density  and  the 
electron  production  rate  at  a  given  height  h.  Case  (a)  would  result  in  the 

lag  time  between  the  growth  of  the  production  rate  and  the  increase  of 
the  electron  density  to  its  equilibrium  values;  Case  (b)  would  result  in  the 
nwytmiyii  lag  time  between  the  reduction  of  the  production  rate  and  the  consequent 
decay  of  the  electron  density. 

Case  (a); 


q(t)  =.  0,  N(h)  =0;  f  <  0; 
q(t)  =  q(h),  N(h)  -  N(h,t);  f  «  0. 
m  this  case  the  production  rate  is  zero  prior  to  a  certain  reference  time, 
t  s  0.  For  times  t  2  0  the  production  rate  proceeds  at  its  full  Intensity  as  a 
function  of  height  as  it  would  be  computed  by  the  method  given  earlier.  N(h) 
Increases  from  zero  and  in  some  characteristic  time,  which  is  a  function  of 
height,  assumes  a  value  arbitrarily  close  to  its  theoretical  equilibrium  value. 
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(^.2) 


Case  (T?); 

q(t)  =  q(h),  N(h)  =  N^Ch);  t  <  0; 

q(t)  -  0,  N(h)  »  N(h,t);  t  S  0. 

Li  this  case  the  equUlbrlua  condition  Is  assumed  to  hold  between  q  and  N. 

Nq  Is  the  equilibrium  value  of  the  electron  density  corresponding  to  q  at  a  given 
height,  h.  At  t  i  0,  q  >  0,  the  primary  production  of  electrons  suddenly  ceases 
and  N  decays  from  Its  equilibrium  value. 

The  eqviatlons  to  be  solved  subject  to  the  above  initial  conditions,  Case  (a) 
and  Case  (b),  are  (3*3)  and  (3*5), 

dN/dt  =  q  -  -  Tjn(0g)N  +  anN"  +  pSN",  (3-3) 

dSf/dt  =  nn(02)N  -  unN"  -  -  pSN",  (3-5) 


with  N  +  N  «  N  . 


If  X,  the  negative- Ion-electron  ratio,  could  be  considered  a  constant, 
special  solutions  of  simple  analytical  form  exist  for  the  above  system  of  equa¬ 
tions  and  the  Initial  conditions  postulated.  However,  an  assumption  such  as  this 
Is  not  Justified  without  some  evidence  as  to  the  range  of  heights  over  which  it 
Is  valid.  If  Indeed  It  Is  valid  at  all.  Fortunately,  some  evidence  bearing  on 
the  behavior  of  X  can  be  obtained  from  Eqs.  (3*3)  and  (3*^)  under  the  initial 
conditions  of  Case  (b).  Considering  the  special  conditions  of  Case  (b)  first, 
one  can  write  the  system  (3*3)  and  (3 •5)  as  follows: 

dN/dt  =  -OjjNN^  -  nn(0g)N  +  «nN"  +  pSN"  ;  (U.3) 

dN'/dt  »  nn(0g)N  -  anN"  -  o^nV  -  pSN"  ;  (3 *5) 

N”  ■  XN  and  n'*’  =  (1  +  X)N  for  t  >  0,  and  for  t  «  0 

X(t  .  0)  -  X^  .  i|n(0g)/(itn  +  pS)  .  (3.11) 

Multiplying  (^.3)  by  n',  (3*3)  by  N,  subtracting  the  equations  and  dividing  the 
difference  by  N^,  one  obtains  the  following  for  the  left-hand  side  of  the  result¬ 


ing  equation: 
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N(dN"/dt)  -  N"(dN/dt) 

1? 


The  complete  equation  hecomes: 

dX/dt  =  (Op  -  +  Tjn(02)(l  +  X)  -  (tm  +  pS)X(l  +  X)  .  (**.4) 

If  the  first  teim  on  the  rl^t-hand  side  of  Eq.  (4.14'),  i.e.,  [(dtp  -  aj^)Xl("'’] 
is  neglected^  which  is  an  approximation  similar  to  that  made  by  Bailey  in  his 
determination  of  the  equilibrium  value  of  X,  the  equation  can  be  written  in  the 
following  form: 


dX/dt  »  iin(02)(l  +  X)(l  -  X/Xg)  ,  (4.5) 

for  t  >  0,  with  the  initial  condition  X  =  Xg  »  +  pS)  at  t  »  0.  One 

can  see  innnedlately  that  the  differential  equation  (4.5)  with  this  initial 
condition  on  X  ccm  only  be  satisfied  by  the  trivial  solution 

X  =  Xg  ,  for  t  6  0, 

for  if  X  decreased  from  its  equilibrium  value,  the  right-hand  side  of  Bq.  (4.5) 
would  be  positive .  This  involves  a  contradiction,  as  a  decrease  In  X  with  time 
demands  a  negative  first  derivative  which  is  inconsistent  with  the  differential 
equation  (4.5).  The  converse  (X  Increasing  from  Xg  with  time)  also  involves  the 
same  contradiction.  To  an  approximation  of  the  order  of  that  made  by  Bailey,  X 
is  constant  in  tine  over  heifidit  raxiges  of  Interest  for  a  decaying  electron 
density. 

Uhfortiinately,  a  single  analytic  solution  corresponding  to  the  initial 
conditions  of  Case  (a)  cannot  be  obtained  due  to  the  Inhomogeneity  introduced 
into  the  differential  equation  for  X  by  q.  However,  the  solution  corresponding 
to  Case  (b)  indicates  that  X  changes  little  for  large  changes  in  the  electron 
production  rate  and  the  electron  density.  It  would  seem  reasonable  to  asexane 
that  X  can  be  considered  a  slowly  changing  parameter  for  an  increasing  electron 
density  as  well  as  a  decaying  one. 
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To  obtain  a  differential  eqxiatlon  for  the  electron  density,  one  adds  (3 >3) 
and  (3*^)>  remeidberlng  that  «  XN, 

(d/dt)  (1  +  X)N  =  q  -  (1  +  X)(ojj  +  Xa^)H^  .  (4.6) 

Because  X  Is  considered  a  constant  and  equal  approximately  to  X^, 

dN/dt  «  q/(l  +  X)  -  (Ojj  +  Xa^)lJ^  .  (4.7) 

under  the  conditions  assximed  for  Case  (b),  q  »  0  for  t  k  0.  Direct  integration 
of  Eq.  (4.7)  with  q  *  0  yields: 

N  -  N^Cl  +  NQ(ap  +  Xaj)t]  .  (4.8) 

The  electron  density  decreases  essentially  as  t"^. 

To  obtain  a  measure  of  the  rate  of  recoii&>lnation  as  a  function  of  height, 
one  may  inquire  as  to  the  time  required  for  the  electron  density  to  decay  to  l/lO 
of  Its  equilibrium  value. 


(4.9) 


O.Ul^  -  Jijil  +  N^(ajj  +  Xa^)Tl  , 

Tj^Ch)  =  9/N^(h)(ajj  +  Xa^)  , 
where  Tj^(h)  will  be  termed  the  “decay  time."  Tj^(h)  vs  h  is  shown  In  Pig.  12. 

A  more  interesting  quantity  than  the  decay  time  Is  the  rise  time.  For 
the  conditions  previously  set  forth  in  Case  (a),  i.e.,  q  »  0,  N  =  0  for  t  <  0; 
q  s  q(h),  N  >  n(t)  for  t  S  0,  what  period  of  time  must  elapse  for  N(t)  to 
Increase  from  zero  to  9/l0)  of  its  equilibrium  value?  As  mentioned  previously, 
we  have  no  anedysls  in  X  guaranteeing  Its  approximate  constancy  as  we  have  for 
Case  (b).  However,  such  an  approximation  seems  reewonable  in  light  of  the 
behavior  of  X  for  relatively  large  changes  in  If  and  q  idilch  occur  in  the  case  of 
decay  from  equilibrium.  Accordingly  one  Integrates  Eq.  (4.7)  with  q  »  q(h). 

O  O 

The  two  identities,  (l)  sech  x  «  1  -  tanh  x  axid  (ll)  (d/dx)  tanh  x  » 

O 

sech  X,  suggest  a  particular  solution  of  the  form  N  >■  A  tanh  (Bt).  inserting 
this  into  Eq.  (4.7)  and  noting  Identity  (l)  one  sees  that  the  conditions 

((1  +  X)(ap  +  Xa^)/qlA^  =  1  , 


and 
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(1  +  X)/q  AB  =  1 


nust  hold. 

Solving  these  simultaneously^  one  acquires  a  particular  solution  for  Eq. 


N  =  [q/(l  +  tanh  [<l(ajj  +  Xa^)/(l  +  X)]^^t,  (4.10) 

vhlch  satisfies  the  initial  condition,  N  °  0  for  t  =  0,  and 

K  -*  [q/(l  +  ^)(a]j  +  Xofj^)]^^  =  Nq  i  for  t  . 

3be  rise  time  vlll  be  defined  as  the  time  required  for  N  to  attain  a  (9/l0)E  , 


T2(H) 


L 


(1  +  X) 


q(h)(ap  +  Xa^) 


(1  +  X)  nl/2 


qCOjj  +  Xa^) 


-  J/  ii 

$ 


(4.11) 


vhere  q  and  X  (eqxiillbrlxim  value)  are  functions  of  h.  Tg(h)  is  the  "rise  time" 
and  as  a  function  of  h  is  shown  in  Fig.  12. 

It  vlll  be  noticed  that  both  the  decay  and  rise  times  are  in  the  order  of 
20-4o  minutes  at  heights  vhlch  are  Important  for  the  attenuation  of  cosmic  radio 
noise.  The  duration  of  the  steady  particle  bombardment  vhlch  can  be  considered 
a  characteristic  time  for  the  phenomena  is  measured  in  days.  Moreover,  the  catses 
conslc  red  here  are  not  duplicated  in  nature—they  are  idealized  extremes.  Con¬ 
sequently  it  would  seem  that  the  quasi-equilibrium  approximation  is  a  suitable 
one  in  the  present  analysis. 


CHAPTER  V 


Cosmic  noise  Absorption 

From  the  phenomenological  view  point,  the  effect  of  the  solar  pazi^lcle 
bombardment  on  the  atmosphere  is  the  creation  of  a  finite  conductivity  in  a 
region  which,  under  ordinary  circioBstances,  behaves  as  a  pure  dielectric. 

Noimally  the  non- ionized  region  belov  the  D  layer  eidilblts  an  index  of  refraction 
slightly  greater  than  unity  vhlch  veurles  directly  as  the  air  pressure  emd  is 
frequency- independent  over  a  very  large  frequency  range  (essentially  from  zero 
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to  10^^  cycles/aec).  The  presence  of  a  concentration  of  free  electrons  serves 
to  depress  the  index  of  refraction  to  veLLues  less  than  unity  and  to  Introduce 
a  dissipative  process  into  the  propegation  of  radio  waves  through  the  region. 
Purthemore,  hath  the  index  of  refraction  and  the  absorption  are  frequency- 
dependent.  The  region  in  its  ionized  state>  in  contrast  to  its  normal  state,  is 
dispersive.  For  waves  of  frequencies  in  the  neighborhood  of  and  less  than  the 
electron  gyrofrequency  (~  1.5  mc/sec),  the  anisotropy  introduced  into  the 
ionized  atmosphere  by  the  presence  of  the  terrestrial  magnetic  field  can  have 
important  effects  on  their  propagation.  At  these  frequencies  the  medium  is 
birefringent  and  the  amount  of  absorption  depends  on  the  mode  of  propagation. 

It  is  necessary,  therefore,  to  examine  the  severed  effects  that  the  solar  penrbicle 
bombardment  may  possibly  induce  in  the  propagation  of  radio  waves  in  the  height 
range  of  interest.  The  relative  infportance  of  the  various  propagation  effects 
possible  depend  primarily  on  the  radio  frequency  under  consideration  and  the 
magnitude  of  the  free-electron  density  and  the  electron-neutral  collision  fre¬ 
quency. 

Ibe  radio  frequency  chosen  for  the  analysis  presented  in  this  thesis  is 
27.6  mc/sec.  Ibis  was  the  operating  frequency  of  the  riometers  situated  at 
various  locations  in  the  earth's  polar  latitudes  to  study  the  enhancement  of 
the  absorptive  power  of  the  atmosphere  following  various  solar  and  geophysical 
events  (25).  This  choice  of  frequency  has  the  obvious  virtue  of  allowing  a 
direct  comparison  to  be  made  between  calculated  results  and  available  rlometer 
data.  It  also  permits  a  considerable  simplification  in  the  calculation  of  the 
specific  absorption  of  the  atmospheric  region  principally  affected  by  the  solar 
particle  bonbardment. 

The  earth's  magnetic  field  at  Farewell,  Alaska,  the  location  of  lowest 
geomagnetic  latitude  considered  in  this  study,  is  inclined  from  the  zenith  at 
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an  angle  of  l6  degrees.  At  all  other  locations  under  consideration  the  inclina¬ 
tion  is  leas  than  this.  Consequently  the  incoming  cosmic  noise  is  being  propa¬ 
gated  quasi- longitudinally  with  respect  to  the  earth's  field.  For  this  type  of 
propagation  at  27*6  mc/sec,  magneto- ionic  double  refraction  can  be  neglected. 

Both  inodes  of  propagation  follow  similar  paths  and  negligible  difference  exists 
in  their  respective  absorptions. 

Peak  electron  densities  produced  by  solar  particle  bontbardments  are,  at 
the  most,  of  the  order  of  that  of  the  D  layer.  The  index  of  refraction  for  radio 
waves  having  frequencies  in  the  neighboxhood  of  30  mc/sec  in  this  region  is 
approximately  equal  to  one.  As  a  result  of  this,  only  the  elementary  theory  of 
radio-wave  absorption  due  to  dissipative  collisions  between  electrons  and  neutrals 
in  a  all^tly  ionized  plasma  is  required  to  stucty  the  effect  of  the  solar  pcutlcle 
boadbardaent  on  the  upper  atmosphere.  The  absorption  caused  thereby  is  strictly 
nondevlatlve. 

The  specific  absorption  in  decibels  per  km  is  given  by 
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- - X  10  ,  the  square  of  the  plasma  frequency, 
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electron  density  (electrons/cm^), 
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electronic  charge  >  1.6  X  10  coulon^s, 
electron  mass  >  9*11  ^  10  km, 

permittivity  of  free  space  »  8.85  x  10~^  farads/meter, 

Q 

velocity  of  light  >  3  x  10  io/b^c, 
index  of  refraction  1, 
electron-neutral  collision  frequency, 

angular  frequency  of  the  radio  waves  -  I.63  x  10  radians/sec. 
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For  coniputatlonal  purposes  the  specific  absorption  can  be  written: 

X(h)  -  k.6  X  lo\N(h)v(h)/[cD^  +  v(h)^])  db/jon  .  (5-2) 

N(h)  is  the  equilibrium  electron  density  in  the  atmosphere  height  range  30-90 
km  obtained  in  the  manner  described  in  the  preceding  sections.  The  values  used 
for  the  electron  collision  frequency  in  this  thesis  were  those  given  by  Nicolot 
(9).  Nlcolet  arrived  at  the  values  presented  in  his  paper  by  deducing  the  effec¬ 
tive  scattering  cross  section  for  electrons  in  air  from  data  obtained  esqperlmen- 
tally  by  Kane  (26).  Kane,  using  a  rocket-borne  transmitter^  measvired  the  differ¬ 
ence  in  absorption  between  a  7*75  mc/sec  signal  and  its  6th  bamonlc  traversing 
the  atmosphere  durltig  a  polar  blackout.  The  launching  was  made  at  Ft.  Churchill^ 
Canada.  The  data  acquired  by  Kane  is  presented  in  his  paper  as  a  collision  fre¬ 
quency  profile  for  the  atmospheric  height  range  60-80  km.  Nlcolet  made  a 
theoretical  extension  of  this  range  to  include  the  30-90  km  range  needed  for  the 
detexmlnatloo  of  the  total  vertical  absorption  due  to  the  particle  bombardment. 
The  values  of  the  collision  frequency  deduced  by  Nlcolet  from  Kane's  more  limited 
results  are  less>  by  a  factor  of  three>  than  those  values  obtained  from  the  col¬ 
lision  cross  section  of  slow  electrons  with  nitrogen  determined  from  a  laboratory 
ergperlment  of  Anderson  and  Goldstein  (27).  The  atmosphere  collision  frequency 
profile  obtained  from  Anderson's  and  Goldstein's  laboratory  measurement  was  used 
by  Bailey  in  his  paper  on  solar  cosmic  ray  Induced  absorption. 

The  results  of  the  present  calcnilatlon  are  given  in  Fig.  10  which  gives  the 

specific  attenuation  as  a  function  of  height  for  deytlae  and  nighttime  condl- 

tlcos  due  to  a  flux  of  protons  (j^(ll2  Mev)  «  100;  7^  «  $],  and  a  heterogeneous 

flux  of  protons  [j  (112  Her)  »  100;  7^  >  and  a  particles  [J^(119.^  Mev)  »  10; 
p  P  « 

7q  »  ^.73]  of  nua^r  ratio  of  10  to  1.  Contours  of  the  specific  attentuation  as 
a  function  of  geomagnetic  cutoff  due  to  protons  under  daytime  conditions  for 
constant  heights  are  given  in  Figs.  20,  21,  and  22. 

The  total  vertical  attenuation  for  a  given  geomagnetic  particle  cutoff 
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energy  Is  obtained  from  the  integral  of  the  specific  attenuation  over  the  height 
range  of  the  region  of  the  atmosphere  rendered  absorbing  by  the  solar  particle 
Influx 


.  oh  N(h)v(h) 

A(db)  »  4.6  X  10  /  -s - 5  dh  . 

'^h^  to  +  v(h) 


(5-3) 


The  limits  h^^  and  hg  have  been  taken  to  be  30  Ion  and  90  km  respectively.  For 
heights  greater  than  90  km  the  particle  energy  loss  in  the  atmosphere^  due  to 
Its  very  low  density,  Is  negligible  while  primary  electron  production  due  to 
solar  electromagnetic  radiation  becomes  Important.  For  heights  less  than  30  km, 
recoinblnation  is  so  rapid  that  an  ambient  electron  density  is  vlrtuedly  nonexis¬ 
tent. 


Ihe  relative  Increase  in  vertical  absorption  of  cosmic  noise  at  2J.6  mc/sec 
during  a  solar  proton  event  over  that  due  to  the  atmosphere  under  nozmal  condi¬ 
tions  was  calculated  for  Stormer  proton  cutoff  energies  of  261.3,  112,  54,  and 
20.6  Mev  corresponding  to  geomagnetic  latitudes  6l.6,  64.7,  66. 9>  and  68.6 
degrees  resi)ectlvely.  ^e  absorption  as  a  function  of  proton  flux  is  shown  in 
Fig.  11  for  power- lav  proton  spectra  with  7  -  3#  5/  and  J  for  day  and  night. 

The  variations  of  the  absorption  with  the  proton  Stormer  cutoff  energy  (E^)  for 
various  total  flux  values  normalized  to  >  112  Mev  (College,  Alaska)  are  given 
for  daytime  and  nl^ttlme  conditions  for  the  same  three  values  of  7  as  above. 
These  results  are  presented  in  Fig.  23  throui^  26.  Finally,  the  absorption  is 
shown  as  a  function  of  geomagnetic  latitude,  day  and  nlgbt  in  Fig.  29. 


CHAPTER  VI 


The  Effect  .of  the  Earth’s  Field  on  the  Proton  Stopping  Heights 

The  earth's  magnetic  field  exercises  an  influence  of  central  Importance  in 
the  establishment  of  the  electron  production  rate  hel£^t  profile  at  a  given 
geomagnetic  latitude.  Moreover,  in  its  role  as  a  rigidity  spectrometer,  the 
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terrestrial  field,  In  conjunction  vlth  the  molar  particle  spectrum,  determines 
the  latitude  variation  of  the  cosmic  radio  noise  absorption  over  the  polar  cap. 
Yet,  In  the  evaluation  of  the  integral  for  the  electron  production  rate,  Eq. 

(2.9)>  the  earth's  field  Is  assumed  to  exert  a  minor  influence  on  the  particle 
trajectories  In  air— in  fact,  an  effect  svifficlently  small  to  be  negligible. 
Accordingly,  the  solar  particles  could  be  thoxight  of  as  traversing  rectilinear 
paths  through  the  atmosphere,  as  was  assiuned. 

At  first  glance  It  seems  paradoxical  that  the  Influence  of  the  terrestrial 
field  is  of  almost  paramount  Importance  in  one  stage  of  the  analysis  and  can  be 
neglected  In  the  other.  The  reason,  however,  resides  in  the  scale  of  the 
phenomena  Involved  In  the  two  stages  of  the  calculation.  The  geomagnetic  lati¬ 
tude  variation  of  the  cosmic  z«.y  cutoff  rigidity  derives  from  the  over-all 
dipole  structure  of  the  tenrestrlal  magnetic  field  and  Its  influence  on  a  pcu> 
tide  over  a  relatively  long  path  through  space.  On  the  other  hand,  the  Influence 
exerted  by  the  earth's  field  on  a  particle  slowing  down  In  the  atmosphere  depends 
on  the  local  structure  of  the  field  end  Its  interaction  with  the  peu^icle  over  a 
relatively  short  path  through  the  atmosphere.  Thus,  while  considered  in  the 
large,  the  effect  of  the  terrestrial  field  on  the  over-all  peu^lcle  trajectory 
might  be  great.  It  can  also  be  true  that  under  favorable  conditions  the  influence 
of  the  field  on  the  terminal  portion  of  the  particle  trajectory  can  be  entirely 
neglected. 

The  local  terrestrial  field  in  the  north  polar  regioxis  Is  directed  downwfuxl 
and  Is  inclined  from  the  ecash's  normal  at  a  rather  small  emgle.  It  can  be 
considered  to  be  of  uniform  intensity  over  the  height  range  30  to  90  km  above 
the  earth's  surface  at  the  geographic  locations  for  which  the  calculations  con¬ 
tained  In  this  thesis  were  performed.  Presented  In  the  following  table  are  the 
local  magnetic  field  strength,  the  field's  inclination  from  the  zenith  direction. 
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and  the  mlnlioiun  allowable  particle  energy  (cutoff)  for  each  of  the  geographic 
locations  (26)  of  Interest  In  this  thesis. 


Location 

Local  Field  Strength 
(gauss) 

Magnetic  Inclination 
(degrees) 

Cutoff  Energy 
(Mev) 

Pt.  Barrow 

0.57 

9.6 

20.8 

Ft.  Yukon 

0.57 

11.7 

54 

College 

0.55 

12.0 

112 

Farewell 

0.55 

15.8 

261.5 

It  will  be  recalled  that  the  quantity  Y(E;d,(p^h),  appearing  In  the  Integreil 

-2  -1  -1 

of  Sq.  (2.2) t  which  represents  the  number  of  peui^lcles  cm  ster  sec  having 
energies  In  d£  at  E  at  height  h,  was  obtained  from  the  dlfferentleil  energy  spec¬ 
trum  of  the  aolar  particles  Incident  on  the  top  of  the  atmosphere. 

Y(E;0,9>b)  «  N(B^)  .  (2.3) 

The  connection  between  E  at  height  h  and  E^,  the  Initial  energy^  was  obtained 
through  the  known  range-energy  relationship  (U)  for  protons  In  air.  This  method 
requires  a  knowledge  of  the  particle  trajectories ,  and.  In  the  absence  of  a 
loceQ.  magnetic  field,  the  trajectories  are  known  to  be  straight  lines.  If  the 
magnetic  field  were  absent,  a  vertically  incident  particle  would  have  a  range  In 
air  numerically  equal  to  the  atmospheric  pressure  at  Its  stopping  height.  A 
particle  of  oblique  Incidence,  traveling  In  a  direction  Inclined  at  an  angle  6 
from  the  zenith  and  havelng  a  remge  In  air,  R,  would  stop  at  a  vertical  hel^dit 
in  the  atmosphere  corresponding  to  the  pressure  R  cos  0.  Tills  Is  not  true  In  the 
presence  of  a  magnetic  field.  The  field  Introduces  a  curvature  into  the  particle 
trajectory  the  degree  of  which  depends  on  the  magnitude  of  the  particle  mooientimi, 
the  strength  of  the  field  and  direction  of  the  magnetic  field  with  respect  to  the 
particle  momentum  vector.  The  more  energetic  the  particle,  the  more  rigid  Is  Its 
trajectory  and  the  more  closely  Is  the  field- free  situation  approximated.  The 
values  given  In  the  previous  table  for  the  quantities  of  Importance  In  this  res- 
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pect  seem  to  Indicate  that  the  Influence  of  the  local  field  could  be  ignored 
safely.  However,  this  assumption,  if  unwarranted,  would  introduce  an  error  in 
the  electron  production  rate  calculation  that  could  invalidate  the  results  of 
the  entire  study.  For  this  reason  a  separate  Investigation  was  made  of  the 
behavior  of  a  set  of  particles  slowing  do\m  in  the  atmosphere  vmder  the  influence 
of  the  loccd  magnetization  field  characteristic  of  the  four  geographic  locations 
of  Interest  and  having  Initial  conditions  which  would  Insure  that  their  actual 
trajectories  would  exhibit  the  maximum  deviation  from  rectlllnearlty.  13ie 
results  of  this  investigation  are  presented  in  Fig.  13 •  Four  groups  of  pcurtlcles 
were  considered,  each  group  initially  havlxig  the  minimum  allowed  energy  at  one  of 
the  geographic  locations  for  which  the  calculation  was  made.  The  particles  com¬ 
prising  each  group  were  divided  into  subgroups,  each  subgroup  being  assigned  an 
initial  angle  of  incidence  with  respect  to  the  zenith.  Because  the  local  mag¬ 
netic  field  through  Its  own  inclination  from  the  zenith  destroys  the  rotational 
syimetry  of  the  problem,  the  particles  comprising  a  subgroup  were  assigned 
various  initial  azimuthal  directions.  The  straight  horizontal  lines  in  Fig.  13 
represent  the  stopping  heights  of  particles  having  initial  energies  indicated  in 
the  figure  and  incident  at  the  given  zenith  angles  under  the  condition  of  zero 
magnetic  field.  The  sinusoidal- like  curves  represent  the  stopping  heights  of 
the  same  particles  as  a  function  of  their  initial  azimuthal  directions  at  a 
given  zenith  angle  of  incidence  under  the  Influence  of  the  eeurth's  field  associ¬ 
ated  with  the  location  indicated.  For  Farewell,  Ft.  Yukon,  and  Ft.  Barrow  the 
stopping  height  cvurve  for  each  zenith  angle  subgroup  was  drawn  through  points 
representing  the  stopping  heights  of  four  particles  having  initial  azimuth 
directions  of  0,  90,  l8o,  and  270  degrees  respectively.  The  plane  of  zero 
azimuth  is  defined  as  that  containing  the  zenith  direction  and  a  unit  vector 
pointing  outward  but  parallel  to  the  earth's  field.  Trajectories  were  calculated 
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for  particles  of  eight  Initial  azimith  directions  for  each  zenith  angle  for 
College  to  assure  the  smoothness  of  the  stopping  height  curves  for  that  location. 
The  equations  of  motion  of  the  particles  and  the  mathematical  details  of  this 
study  are  given  In  the  Appendix. 

For  each  of  the  four  geographic  locations  considered,  the  stopping  heights 
of  those  particles  having  the  m'<n<wmn  allowed  energy  at  each  location  were  com¬ 
puted.  In  all  cases  the  assumption  of  rectlllnearlty,  carried  to  the  limit  of 
particles  stopping,  appears  to  break  down  at  large  zenith  angles,  but  holds  up 
quite  well  out  to  60  degrees.  As  might  be  expected.  Farewell  represents  the 
worst  case,  as  here  the  magnetic  field  has  a  substantial  horizontal  component, 
being  Inclined  16  degrees  from  the  zenith.  It  must  be  borne  In  mind  that  the 
ln«^v^lllllm  effect  produced  by  solar  protons  on  the  atmosphere  with  respect  to 
IncreeuBlng  Its  opacity  to  cosmic  radio  noise  occurs  at  heights  considerably  In 
excess  of  the  particle  stopping  heights,  with  the  exception  of  Pt.  Barrow.  This 
latter  location,  however,  exhibits  trajectories  having  the  least  deviation  from 
the  stopplxig  heights  obtained  from  the  field- free  analysis. 

This  analysis  cannot  In  itself  provide  a  direct  estimate  of  the  error  in 
the  final  results  Incurred  by  the  assumption  of  rectilinear  pcurtlcle  trajectories 
In  the  atmosphere.  But,  because  It  shows  that  deviations  from  rectlUnesurlty 
are  not  prohibitively  large  for  the  trajectories  of  those  particles  having 
Initial  conditions  and  path  lengths  most  imfavorable  to  this  assumption,  it  does 
Indicate  that  the  over-all  error  introduced  Into  the  final  results  thereby 
cannot  be  large,  and,  in  fact,  ml^t  well  be  entirely  negligible  considering 
the  uncertainties  surrounding  other  aspects  of  the  calculation  of  the  cosmic 
noise  absorption.  In  addition,  it  should  be  noted  that  the  earth's  cosmic  ray 
shadow  establishes  foibldden  sectors  affecting  particles  incident  from  large 
zenith  angles  (l?)  thus  reducing  the  number  of  particles  whose  trajectories 
violate  the  assumption  of  rectiUnearlty. 
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OHAPTER  VII 

DlscuBslOD  of  the  Results 

Ibere  are  tvo  distinct  aspects  from  vhlch  one  can  consider  the  various 
phenomena  attendant  \ipon  a  solar  particle  bombardment  of  the  atmosphere.  From 
the  purely  local  point  of  view  one  may  examine  the  various  Intezaction  processes 
occurring  In  the  atmosphere  as  a  result  of  the  Influx  of  solar  cosmic  rays  as 
veil  as  the  distribution  with  height  of  the  electron  production  rate,  the  elec¬ 
tron  density  and  the  specific  attenuation.  Of  special  Interest  are  the  various 
“layer  heights"  or  heights  of  maximum  value  of  these  quantities.  A  second  and 
equally  Interesting  aspect  of  the  over-all  particle- Induced  absorption  jbenomenon 
exists  on  a  geographical  scale.  Here,  one  considers  the  geomagnetic  latitude 
variation  of  the  total  attenuation  of  the  cosmic  radio  noise,  and  the  geomagnetic 
cutoff-energy  contours  of  the  electron  production  rate,  electron  density,  and 
specific  attenuation  for  a  set  of  given  atmospheric  heights. 

It  can  be  stated  at  the  outset  that  a  distinguishing  characteristic  of  the 
local  height  profiles  of  the  varlovis  quantities  of  Interest  (production  rate, 
electron  density,  and  specific  attenuation)  Is  their  virtual  Independence  of 
the  shape  of  the  solar  particle  energy  spectrum.  The  profiles  seem  to  be 
detezmlned  solely  by  the  atmosphere,  while  the  absolute  magnitudes  assumed  by 
the  quantities  are  determined  by  the  size  of  the  particle  flux  In  the  neighbor¬ 
hood  of  the  local  cutoff  enezgy.  On  the  other  hand,  these  same  quantities,  when 
viewed  through  their  cutoff  energy-constant  hel^t  contours,  show  a  strong 
dependence  on  the  shape  of  the  solar  particle  spectrum.  Consequently,  the 
latitude  variation  of  the  total  attenuation  of  cosmic  raidlo  noise  should  be 
deduclble  from  the  primaxy  solar  cosmic  ray  spectrum  provided  a  known  geomagnetic 
cutoff  mechanism  Is  operative.  This  suggests  the  attractive  possibility  of 
obtedning  the  solar  particle  energy  spectrum  directly  from  Ionospheric  absorption 
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data  gathered  at  various  geoiuagnetlc  latitudes.  Before  this  can  be  done  with 
confidence,  however,  all  details  of  the  absorption  phenomenon  on  the  local  level 
must  be  veil  understood.  This  applies  In  particular  to  the  over-all  reconiblna- 
tlon  process  through  whose  agency  the  ambient  electron  density  Is  established. 

The  numerical  values  of  the  rate  coefficients  of  processes  known  to  be  li^portant 
are  controversial,  and  It  veil  nay  be  true  that  every  Important  process  has  not 
been  accounted  for.  It  Is  to  emphasize  these  differences  that  some  of  the 
results  obtained  by  Bailey  In  his  study  on  cosmic  ray  Induced  absorptloi  (12) 

80%  compared  In  this  thesis  with  those  obtained  from  the  present  study.  More¬ 
over,  the  variation  of  the  low-energy  cutoff  of  the  cosmic  ray  spectrum  with 
latldude  must  be  known.  Does  a  Stormer-type  mechanism  operate  following  a 
solar  flare  before  the  onset  of  the  Inevitable  magnetic  stom,  or  Is  the  geomag¬ 
netic  field  always  grossly  perturbed  during  the  whole  course  of  such  an  event? 

To  obtain  the  solar  peortlcle  spectznmi  from  absorption  data  one  must  be  able  to 
predict  the  latitude  behavior  of  the  low-energy  cutoff  as  well  as  account  for 
the  various  atmospheric  processes  mentioned  above. 

The  Local  Aspect  of  the  Terrestrial  Effects  of  a  Solar  Flare 

We  begin  the  discussion  of  the  results  of  the  analysis  described  In  pre¬ 
vious  sections  of  this  thesis  with  a  consideration  of  a  solar  event  as  viewed 
locally,  with  enqphasls  at  present  on  the  hel^t  profiles  of  the  various  important 
quantities  connected  with  the  event.  The  rate  at  idiich  electrons  are  produced 
in  the  atmosphere  by  a  particle  bombardment  is  a  proper  first  consideration. 

Tbe  curves  representing  the  electron  production  rate  presented  In  this  thesis 
obviously  depend  on  the  magnitude  of  the  bombarding  flux  assumed.  The  particular 
value  chosen  for  noimadlzlng  the  various  solar  proton  energy  spectra  for  which 
calculations  were  done  was  chosen  as  being  representative  of  the  flux  magnitudes 
encountered  in  the  aftermath  of  a  large  proton  event.  The  locale  chosen  for 
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special  study  was  College,  Alaska  (geoaagnetlc  latitude  64.7°;  geonagnetlc 
cutoff  E  >  112  Mev).  J(ll2  Mev)  «  100  protona/em^'sec-ster  agrees  closely  In 
magnitude  with  the  actual  Integral  spectnim  obtained  by  Winkler  et  al.  (15)  at 
Minneapolis  following  the  sudden  ccnneneement  of  a  magnetic  storm  associated 
with  the  polar  cap  event  of  July  l4,  1959*  The  spectrum  obtained  as  reported  by 
these  authors  was 

j(E)  =  1.05  X  10®E"^‘^  protons/cm^-sec-ster,  (7*1) 

f or  88  <  E  <  300  Mev. 

At  College,  Alaska,  the  flux  was 

J(112)  «  (1.05  X  10®)/(112)^*^  »  119.8  protons/cm^-ster-sec.  (7.2) 

The  first  result  of  interest  obtained  from  this  study  has  already  been 
mentioned,  that  Is,  the  seeming  lack  of  dependence  of  the  height  profile  of  the 
electron  production  rate  on  the  cosmic  ray  energy  spectrum.  !Ihe  reason  for  this 
becomes  apparent  \dien  one  glances  at  Fig.  1,  the  energy- loss  curve  for  protons. 
Low-energy  particles  are  more  effective  Ionizing  agents  than  those  of  high 
energy,  and  there  are  more  of  them  If  a  power  law  of  any  negative  order  approxi¬ 
mates  the  particle  energy  spectrum.  Figure  2  displays  this  result.  Ihe  electron 
production  rate  due  to  bombarding  protons  of  three  different  differential  energy 
spectra  are  shown  for  a  geonagnetlc  proton  energy  cutoff  of  112  Mev.  The  curves 
presented  peak  strongly  at  about  32  km,  the  atmospheric  height  corresponding  to 
the  stopping  depth  of  112  Mev  protons  (10.7  gn^cm^).  The  curves  fall  steeply 
with  increasing  helj^t  from  their  peak  valvms  (which  are  almost  equal  due  to 
the  normalization  Imposed  on  the  Integral  spectra)  until  a  hei^t  of  about  58  km 
la  reached.  Thereupon  the  descent,  still  steep,  becomes  a  little  more  gradual. 
This  Is  due  to  the  growing  importance  of  the  laterally  incident  particles.  At 
heights  of  58  km  and  greater,  the  air  Is  so  thin  that  protons  Incident  from  all 
angles  are  contributing  equally  to  the  electron  production  rate. 
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This  Bane  pattern  Is  followed  for  other  values  of  the  lov-energy  cutoff. 

A  strong  peak  is  obtained  at  the  height  corresponding  to  the  stopping  depth  of 
the  lowest  energy  particles.  A  steep  fall-off  follows  the  sane  general  shape  as 
that  given  in  Fig.  2.  This  can  be  seen  in  Fig.  3.  Three  curves  are  displayed 
here:  an  electron  production  rate  profile  due  to  protons  having  a  differential 
energy  spectrum  oST^,  a  profile  due  to  alpha  particles  having  an  Integral 
rigidity  spectrum  Identical  to  that  of  the  protons  but  In  a  peurtlcle  ratio  of 
one  to  ten  with  respect  to  the  protons,  and  a  profile  due  to  the  sum  of  the  two. 
It  will  be  noticed  that  the  alpha-partlcle-electron  production  peak  is  displaced 
from  that  of  the  proton- electron  production  rate  peak  upward  to  30  km.  This 
height  corresponds  to  the  stopping  depth  of  alpha  particles  vertically  incident 
with  an  energy  of  119.3  Mev,  their  geomagnetic  cutoff  corresponding  to  a  proton 
cutoff  of  112  Mev. 

The  alpha-partlclc-electron  production  rate  profile  also  displays  an  Inter¬ 
esting  Interplay  between  the  geomagnetic  cutoff  and  that  due  to  the  atmosphere. 
The  electron  production  rate  represented  by  the  portion  of  Curve  II  (Fig*  3)  to 
the  left  of  Its  peak  Is  entirely  due  to  alpha  particles  Initial  r  having  energies 
In  excess  of  their  geomagnetic  cutoff  value  of  119*3  Ifev.  Because  the  mlnlmun 
energy  that  an  alpha  particle  must  have  to  reach  a  given  depth  for  these  par¬ 
ticles  Is  a  function  of  the  depth  and  zenith  angle  of  Incidence  only,  the  atmos¬ 
phere  cutoff  and  the  particle  spectnim  detennlne  this  part  of  the  electron  pro¬ 
duction  rate  profile.  Hils  seems  to  be  contrary  to  the  statement  made  earlier 
concerning  the  lack  of  dependence  of  the  electron  production  rate  profile  on 
the  particle  spectrum.  It  must  be  remeoibered,  however,  that  this  statement  was 
based  on  the  assumption  of  a  flux  predominantly  made  up  of  protons. 

It  Is  really  more  precise  to  say  that  the  electron  production  rate  height 
profile  Is  Independent  of  the  incident  particle  spectrum  for  that  range  of 
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heights  wherein  Ionization  Is  predoriinently  due  to  peurbicles  that  entered  the 
atmosphere  having  the  geomagnetic  cutoff  energy  characteristic  of  the  given 
latitude.  For  protons  at  College^  Alaska^  these  heights  are  32  km  emd  greater. 
For  heights  less  than  these,  the  electron  production  ray  due  to  protons  falls 
off  sheurply  with  increasing  depth  (decreasing  height)  and  electron  recombination 
is  rapid  and  complete.  Ihe  contribution  to  the  total  cosmic  radio  noise  attenua¬ 
tion  due  to  the  axoblent  electron  density  distribution  that  would  be  spectrum^ 
sensitive  Is  negligible.  Were  ea.pha  peurtlcles  the  primary  ionizing  agent,  there 
would  be  a  local  dependence  of  the  electron  production  rate  on  the  particle 
spectrum. 

The  electron  production  rate,  represented  by  the  portion  of  Curve  II  (Fig.  3) 
to  the  right  of  Its  peak.  Is  due  primarily  to  the  geomagnetic  cutoff  cd.pha 
particles.  Like  that  poi*tlon  of  the  proton-electron  production  rate  profile 
analogous  to  it,  this  part  of  the  production  rate  curve  is  Independent  of  the 
alpha-particle  spectrum. 

The  sum  profile  (Curve  III,  Fig.  3)  shows  em  approximate  twofold  Increase 
In  the  electron  production  rate  resulting  from  the  addition  of  one  alpha  particle 
per  ten  protons  In  the  flux.  This  Is  not  a  drastic  change,  especially  In  view 
of  the  f-act  that  the  ambient  electron  dezuilty  varies  as  the  square  root  of  the 
production  rate  [liq.  (3«13)}>  It  would  seem  that  a  contamination  of  a  proton 
fl\nc  by  alpha  peo^lcles,  unless  it  were  very  heavy,  would  not  be  readily  dls- 
cemable  from  cosmic  radio  noise  measurements . 

The  electron  density  height  profiles  resulting  from  the  action  of  a  solar 
cosmic  ray  bombardment  are  shown  in  Fig.  4.  Features  of  Immediate  Interest  are 
the  locations  of  the  peaks  of  the  curves  presented.  Both  the  daytime  and  nl^t- 
tlrne  curves  attain  their  maxlminn  values  at  considerably  greater  heights  than 
those  at  which  the  maximum  primary  production  of  electrons  occurs.  The  height 
of  peak  electron  density  Is  greater  at  night  than  for  the  day  which  Is  character- 
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Istlc  of  the  befiavlor  of  the  nonnally  occurring  ionosjAiere  layers.  Curve  Il(b) 
(Fig.  4)  represents  the  nighttime  electron  density  profile  due  to  an  alpha-par- 
ticle-proton  flux  and  exhibits  an  irregularity  at  about  52  km  due  to  the 
stopping  aJ.pha  peurtlcles.  Otherwise  the  curves  are  quite  smooth. 

IQie  relative  positions  of  the  daytime  pesdc  and  nighttime  peak  are  opposite 
to  those  given  by  Bailey  (12).  Die  midday  electron  density^  Euscording  to  this 
author^  peaks  at  about  85  km  \dille  the  nighttime  profile  has  its  maximum  at 
about  75  km.  Bailey's  results  and  the  results  derived  from  the  present  analysis 
agree  in  one  respect ^  that  being  the  approximate  constancy  of  the  layer  height 
with  change  in  geomagnetic  latitude.  Oirtner,  Egeland  and  Hulqvlst  {2$),  basing 
their  statement  on  the  results  derived  in  an  earlier  paper  by  Ikilqvist  and 
Ortner  (3^)^  believe  tlact  the  height  of  maximum  electron  density  decreases  with 
deci'easlng  geomagnetic  latitiide.  For  a  geomagnetic  latitude  of  59*7  they  claim 
a  layer  is  formed  by  a  proton  bombardment  having  a  lower  boundeoy  situated  at 
a  height  of  10  km  and  the  "point  of  gravity"  of  the  electron  distribution 
resides  at  15-20  km  above  ground  level.  Ihey  cite  very  low  frequency  and  medium 
freqviency  propagation  observations  to  support  their  position. 

These  differences  with  the  results  of  the  present  study  are  qualitative; 
quantitative  discrepancies  are  even  more  marked.  Of  all  the  limits  in  the 
chain  of  deduction  connection  the  spectrum  of  the  bombarding  solar  particles  to 
the  cosmic  radio  noise  aboszptlon,  the  one  Joining  the  electron  density  profile 
to  the  electron  production  rate  in  the  atmosphere  is  sxirely  the  most  doubtful. 

It  has  been  brought  out  earlier  that  the  electron  recombination  mechanism  in 
the  height  range  30  to  90  km  is  a  complex  one,  and  there  is  disagreement  con¬ 
cerning  the  values  of  the  rate  coefficients  of  the  various  component  processes. 
It  has  also  been  mentioned  that  the  nunerlcal  discrepemcles  appearing  between 
Bailey's  calculation  of  the  negative- ion- to-electron  ratio,  X,  and  those 
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presented  here  arise  from  the  different  rate  coefficients  xised  rather  than  from 
the  differing  methods  employed  in  obtaining  this  quantity.  Hils  is  true  of  all 
the  dlscrepeuicles  displayed  in  Figs.  5,  6,  J,  and  8  vhlch  present  profiles  of 
important  quantities  relevant  to  the  electron  density  distribution.  It  should 
be  pointed  out  that  the  curves  representing  the  effective  recombination  coef¬ 
ficient  as  a  function  of  height  should  be  read  from  right  to  left,  the  abscissa 
being  on  the  top  of  the  figure  for  This  was  done  to  make  Figs.  5  and  6 

more  easily  legible.  On  examination  of  Figs.  3  and  6,  one  observes  em  almost 
constant  order- of -magnitude  difference  in  X,  for  a  given  height,  existing 
between  Bailey's  curve  (B)  and  the  curve  derived  in  the  course  of  the  present 
investigation  (A).  This  discrepancy  reflects  the  order- of -magnitude  difference 
in  the  respective  choices  of  the  values  of  q,  the  attachment  rate  coefficient. 
Whereas  X  as  computed  by  Bailey  is  larger  then  that  derived  in  the  present  work, 
this  is  reversed  in  the  case  of  the  effective  recombination  coefficient. 

Figures  7  and  8  compere  for  day  and  night,  respectively,  the  quantity  (l  + 
which,  in  fact,  if  not  in  name,  is  the  equilibrium  recombixiatlon  coefficient  in 
the  3O-9O  km  height  range.  Ihe  rate  coefficients  used  by  Bailey  dictate  a 
smaller  electron  density  for  a  given  value  of  the  electron  production  rate  than 
possible  for  those  used  in  the  present  study.  This  is  explicitly  shown  in  Fig. 
9,  wherein  the  ratio  of  the  electron  density,  as  would  be  calculated  with 
Bailey's  choice  of  coefficients  (Ng)  to  that  (N^)  computed  with  those  rate 
values  employed  in  the  present  study,  is  displayed.  Except  for  helots  in 
excess  of  79  km  under  daytime  conditions  less  than  unity.  Between  90 

and  60  km,  vhlch  is  the  most  important  height  range  as  far  as  the  actual  absorp¬ 
tion  of  cosmic  radio  noise  is  concerned  (see  Fig.  10),  the  daytime  curve  repre¬ 
sents  nearly  an  order  of  magnitude  difference. 

Figure  10  gives  the  specific  attenuation  (db/km)  profiles  for  cosmic  radio 


noise  at  27.6  mc/sec  corresponding  to  the  electron  density  profiles  of  Fig.  4. 

A  rather  striking  feature  displayed  by  these  curves  is  the  reversal  of  the 
height  positions  of  the  daytime  and  nighttime  peaks  with  respect  to  their  rela¬ 
tive  positions  in  the  electron  density  profiles.  Kie  positions  of  the  peak 
values  of  the  specific  attenuation  occ\ir  at  considerably  greater  heights  than 
given  by  Hulq.vl8t  and  Ortner  (3^)  but  are  slightly  lower  than  the  maxima  of 
Bailey's  curves.  It  is  also  noteworthy  that  the  peak  values  of  the  specific 
attenuation  profiles  occur  at  lesser  hei£^ts  than  those  of  the  corresponding 
electron  density. 

!nie  total  vertical  attenuation  as  a  function  of  the  proton  flux  at  geomag¬ 
netic  cutoff  of  112  Mev  is  given  in  Fig.  11.  These  curves  show  the  almost  neg¬ 
ligible  dependence  of  the  attenuation  on  the  shape  of  the  proton  Integral  spec¬ 
trum.  The  daytime  curves  are  straight  lines  on  a  log- log  scale.  The  slif^t 
forward  bend  in  the  nl£^ttlme  curves  is  due  to  the  dependence  of  X,  the  negattve- 
lon-to-electron  ratio,  on  the  electron  production  rate  for  large  flux  values. 

This  would  not  be  observed  were  X  obtained  by  Bailey's  method.  Ihe  effect, 
however,  is  small  which  demonstrates  the  excellence  of  Bailey's  approximation 
for  X  [Eq.  3.11)]. 

The  Absorption  as  a  Function  of  the  Particle  Ckttoff  Energy 

It  has  been  pointed  out  several  times  in  the  coxirse  of  this  thesis  that 
the  low-energy  part  of  the  solar  cosmic  ray  spectrum  is  responsible  for  most  of 
the  ionization  occurring  in  the  upper  atmosphere  during  a  solar  particle  bombard¬ 
ment.  Consequently,  the  value  of  the  low-energy  cutoff  is  of  decisive  Importance 
in  determining  the  amount  of  absorption  suffered  by  cosmic  radio  noise  at  a  given 
geographic  location. 

The  calculation  of  the  cosmic  radio  noise  attenuation  was  originally  made 
for  cutoff  energies  corresponding  to  the  Quenby-Webber  proton  cutoff  rigidities 


of  four  geogra^dilc  locations  at  which  facilities  were  available  for  meetsurlng 
the  cosmic  noise  attenuation  associated  with  a  solar  event.  The  miwimum  energy 
the  solar  peortlcles  must  possess  to  reach  the  terrestrial  atmosphere  was  thovight 
of  as  being  established  In  a  known  way  by  the  earth's  magnetic  field.  Hence^  a 
particle  cutoff  energy  variation  with  geomagnetic  latitude  would,  through  the 
solan:  cosmic  ray  spectrum,  be  reflected  In  a  characteristic  distribution  of  cos¬ 
mic  radio  noise  absolution.  However,  the  validity  of  the  results  displayed  In 
Figs.  14  through  28  which  give  constant  height  contours  of  the  electron  produc¬ 
tion  rate,  electron  density  and  specific  attenuation  as  well  as  the  total  cosmic 
noise  absorption  eis  a  f\uictlon  of  the  solar  proton  cutoff  energy  for  day  and 
ni£^t.  Is  not  affected  by  the  eissumptlon  made  concerning  the  type  of  cutoff 
operative  dvirlng  an  event.  In  a  later  section  the  measured  geograjhic  distribu¬ 
tion  of  cosmic  radio  noise  attenuation  associated  with  the  solar  flare  of  July 
1959 t  1b  discussed. 

Figures  l4,  1^,  and  l6  give  constant  eight  contours  for  the  electron  pro¬ 
duction  rate  versus  cutoff  energy  due  to  protons  distributed  In  energy  according 
to  the  differential  energy  spectra  Indicated  In  the  figures.  These  curves  can 
be  used  to  facilitate  further  electron  production  rate  calculations  using 
“broken”  spectra,  as  mathematically  Is  but  a  lower  limit,  and  other  chcmges 
necessary  for  such  a  calculation  merely  Involve  adjustments  In  multiplicative 
constants.  For  Instance,  If  one  wished  to  ascertain  the  electron  production 
rate  at  a  given  height  to  a  trial  spectrum  of  the  type 

N(E^)  =  a^Ej3  ^  ^ 

and 

N(E^)  =  OtgE^"^  ,  for  E^  <  E^  S  ^ 


the  curves  could  be  used  In  accordance  with  the  relation 


where  the  expression  represented  by  the  difference  of  the  double  Integrals  on 
the  rlg^t-hand  side  of  the  above  symbolic  equation  can  be  obtained  from  the 
curves  presented  In  Fig.  l4  {eay)  with  a  proper  adjustment  of  constemts.  Ihe 
contours  are  of  sufficient  regularity  that  Interpolation  Is  possible,  l^ls 
offers  a  method  of  obtaining  q(h)  with  a  minimum  of  labor  for  spectra  other  them 
those  chosen  for  the  present  work. 

Figures  I7  througli  22  present  constant  hel|^t  contours  for  the  daytime  elec¬ 
tron  density  emd  specific  attenuation  respectively  as  functions  of  the  solar 
proton  cutoff  energy  for  the  proton  differential  energy  spectra  Indicated. 

Because  both  of  these  quantities,  for  a  given  cutoff  energy,  present  profiles 
which  assume  maximum  values  at  atmospheric  heights  of  some  Intermediate  value 
In  the  height  range  of  Interest,  the  contours  are  drawn  In  such  a  majuier  as  to 
avoid  the  Inevitable  confusion  that  would  result  from  a  conventional  display. 
Contours  for  heights  less  than  the  maximum  are  drawn  frem  left  to  rl£^t— those 
for  heights  greater  than  the  maximum  are  drawn  oppositely.  The  possibility  of 
drawing  these  contours,  as  presented,  demonstrates  the  approximate  constemey  of 
the  height  of  the  layer  maxima  with  geomagnetic  latitude,  were  a  Stoxmer 
mechanism  effective  In  establishing  a  low-energy  cutoff.  This  was  mentioned 
earlier  in  this  discussion  In  connection  with  a  seemingly  contrary  opinion  held 
by  Ortner,  Egelemd,  and  Hulqvist  (29)« 

Figures  23,  24,  and  25,  the  daytime  attenuation  versus  cutoff  energy  curves, 
display  one  of  the  most  striking  results  obtained  In  the  study,  the  rectilinear 


nature  of  these  cur/es  when  plotted  on  log- lop  paper.  Hie  nighttime  curves  are 
not  quite  straight.  This  is  due^  in  their  case,  to  the  dependence  of  the 
negative- ion-to-electron  ratio  on  the  electron  production  rate  at  nlj^it.  Hie 
regularity  of  the  daytime  cvixves  suggests  the  possibility  of  a  simple  empirical 
fomula  which  might  give  the  total  vertical  attenuation  provided  the  cosmic  ray 
flux,  the  cutoff  energy  and  7,  the  slope  of  the  differential  energy  spectrum  (on 
a  log- log  plot)  are  known. 

For  the  daytime  case,  Bailey's  expression  for  X,  the  negative- ion-to-electron 
ratio,  [Eq.  (3.II)]  is  valid.  If  X  is  Independent  of  the  electron  production 
rate,  the  square  of  the  total  vertical  attenuation  is  proportional  to  the  flux. 


where  A  is  the  total  vertical  attenuation  in  db  and  J  is  the  Integral  flux 
(protons/cm^-sec-ster)  above  cutoff.  Moreover,  for  a  given  value  of  J,  the 

«.y  !  2 

total  vertical  attenuation.  A,  is  nearly  proportional  to  E  "  , 


-  I  ♦  P(E  ,r) 

A  «  L(7)E^  ^  ® 


(7.5) 


where  E  is  the  low-encrgy  cutoff  in  Mev,  1(7)  is  a  proportionality  factor,  and 
c 

P(E^,7)  Is  a  slowly  varying  parameter. 

These  facts  suggest  an  investigation  of  the  quantity 


log^QlA\’'/j(ll2)]  , 

where  A  is  total  vertical  attenuation  (db)  of  cosmic  radio  noise  corresponding 
to  a  proton  differential  energy  spectrum  of  W(E)  =  cE  ^  and  a  low-energy  cutoff 
of  E^.  J(ll2)  is  the  corresponding  integral  proton  flux  for  E^  »  112  Mev. 

In  order  to  presezit  numerical  results,  it  is  necessary  to  choose  a  "refer¬ 
ence"  cutoff  energy  for  the  Integral  spectnus.  E  »  112  Mev  is  the  nozmal 
proton  geomagnetic  cutoff  at  College,  Alaska.  However,  the  formula  \dilch  will 
be  presented  presently  does  not  depend  on  this  interpretation  of  the  argument 


of  J(E  )  =  J(112). 
c 


The  fol.lowlijg  empirical  relationship  is  obtained  from  the  curves  presented 
in  Figs.  23,  2k,  and  25: 

log^QlAhj /J{112)]  =  2(7  -  1)  +  W(7,E^)  >  (7.6) 

where  M(7,E  )  is  given  in  Fig.  30  for  7  =  3,  3,  and  7,  a  range  in  E  of  20  S  E  i 
c  c  c 

260  Mev,  and  a  range  of  integral  spectra  10  S  J(ll2)  ^  ^00  protons/cm^-sec-ster. 

Equation  (7*6)  and  Fig.  3O,  to  a  very  good  approximation,  embody  the  results 
of  the  calculations  performed  to  obtain  the  total  vertical  daytime  cosmic  noise 
attenuation  due  to  a  known  proton  flux  at  a  frequency  of  27.6  mc/sec.  The  day¬ 
time  curves  in  Figs.  23,  2k,  sind  25  can  be  derived  from  them,  and  through  inter¬ 
polation,  the  cosmic  noise  attenuation  as  a  function  of  the  proton  cutoff  energy 
for  cosmic  ray  spectra  other  than  those  for  which  the  full  calculation  was  made 
can  be  obtained.  Equation  (7. 6), along  with  the  curves  in  Fig.  30,  can  be  thought 
of  as  representing  an  ultimate  condensation  of  the  results  of  the  calculations 
for  the  total  daytime  attenuation  of  cosmic  radio  noise  at  27.6  mc/sec  due  to  a 
solar  proton  bombardment  of  the  atmosphere. 

Figure  29  shows  the  theoretical  geomagnetic  latitude  distribution  of  cosmic 
radio  noise  attenuation  under  the  assusiption  that  the  terrestrial  magnetic  field 
is  undisturbed  and  that  a  power- law  spectrum  holds  down  to  cutoff.  The  low- 
energy  cutoff  is  obtained  from  the  Quenby-Webber  cutoff  rigidity  formulae  which 
represent  a  further  refinement  of  the  simple  dipole  theory  (16).  The  latitude 
variation  displayed  in  Fig.  29  is  quite  steep.  It  will  be  seen  in  a  following 
section  that  the  measvired  cosmic  radio  noise  attenuation  does  not  display  the 
strong  latitude  variation  shown  in  Fig.  29,  but  instead  tends  to  be  almost 
constant  over  the  whole  polar  cap. 
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CJIAPTER  VIII 


The  Rlometer 

Die  cosmic  radio  noise  attenuation  data^  to  which  the  results  presented  In 
this  thesis  are  to  be  compered,  are  obtained  from  rlometer  measurements  made 
during  a  solar  particle  event.  However,  before  comparisons  are  drawn,  a  brief 
accovuat  should  be  given  of  the  relevant  operational  characteristics  of  the 
recording  Instrument.  For  a  more  complete  description  of  the  rlometer  (Relative 
Ionosphere  Opacity  Meter)  one  Is  referred  to  the  paper  devoted  to  this  instru¬ 
ment  by  Little  and  Lelnbach  (3l)*  For  the  purposes  of  this  thesis  the  following 
data  are  sufficient: 

Receiver  Characteristics 

Center  Frequency  27.6  mc/sec 

Width  of  Frequency  Sweep  100  kc 

R.  F.  Bandwidth  0.3  me 

In  operation,  a  6-kc  exploring  band  Is  swept  at  about  2.5  kc/sec.  Ihls 
band  In  turn  traverses  a  100-kc  search  band  every  40  sec ends.  The  minimum  signal 
strength  recorded  Is  attributed  to  cosmic  radio  noise.  Ihe  sweeping  technique 
described  Is  designed  to  discriminate  against  man-made  narrow  band  interference 
at  the  operating  frequency  of  the  system. 


Antenna- - 3- element  Yagl  with  I/2  power  beamuldth 


{ 


Scf*  E  plane 
110°  H  plane 


DynAwie  Ranf|*-.is  db  (fixed  gain);  30  db  (variable  gain). 

It  is  appeurent  immediately  that  the  single-frequency  cmalysls  performed  for 
this  thesis  is  adequate.  The  spectral  congponents  In  the  operating  band  of  the 
system,  as  given  above,  would  show  entirely  negligible  differential  attenuations. 
The  same  statement  cannot  be  made  with  regard  to  the  total  vertical  attenuation 
cadculated  and  the  attenuation  actually  recorded.  The  relatively  large  beamr 
width  of  the  antenna  employed  pezmlts  the  reception  of  a  significant  amount  of 
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radiation  that  has  traversed  the  absorbing  region  of  the  upper  atmosphere 
obliquely.  Ibis  should  result  in  a  greater  measured  absorption  for  a  given 
particle  flux  than  that  calculated.  !Ibe  difference  can  be  accoimted  for  in  an 
approximate  vay  by  an  integration  over  an  hypothetical  receiving  pattern  that 
simulates  the  actual  pattern  of  the  Yagl.  The  following  brief  analysis  is 
admittedly  rough,  but  demonstrates  the  fact  that  the  results  derived  for  verti¬ 
cal  attenuation  do  not  differ  very  significantly  from  those  measured,  provided 
the  difference  is  due  only  to  the  receiving  pattern  of  the  rlometer  antenna. 

Let  the  noise  power  received  vertically  during  an  undisturbed  period  be  P^, 
and  that  similarly  received  during  flare  conditions  be  P.  Then  their  ratio  is 

P/P^  =  e"®  ,  (8.1) 

where  a  is  the  Integral  of  the  specific  attenuation  over  the  vertical  height 


(in  km). 


lo3 


a  = 


h-(km)  2,. 
2  oj  V 

JL 


dh  . 


(8.2) 


ij^(km)  ttP  + 

Ibe  quantities  in  the  right-hand  side  of  the  above  equation  have  been  defined 
previously.  The  total  relative  vertical  absorption.  A,  is  given  by 

A  »  a(lO  logj^Qe)  =  4.32a  (db)  .  (8.3) 

Oh  the  other  hand,  the  total  relative  absorption,  taking  the  receiving  pattern 
into  account,  is 


r  /  ^  ^  0(e,9)e"*  ®  sin  0  dfi  <kp  -I 

A  -  10  log-oe  -2 - — -75 -  (db)  , 

^  /o  /q  sin  0  d0  dp  -1  (Q 

where  0  is  the  absorption  cross  section  of  the  antenna.  A  \inidlrectlonal  power 
pattern  of  cos^0  approximates  fairly  well  the  Yagl  pattern  with  the  I/2  power 
beamwldths  given  above.  It  has  a  directivity  (gain)  of  4  db  over  isotropic  and 
a  1/2  power  beamwldth  of  9^°  (about  midway  between  the  E  and  H  plane  beamwldths 
of  the  Yagl,  as  quoted  from  Little  and  Lelnbach).  The  shape  of  the  unldirec- 
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tlonal  coslne-squared  pattern  ±b  such  as  to  slightly  enhance  the  effect  of 
obliquity  over  that  due  to  the  Yagl  pattern*  This  can  be  seen  by  an  inspection 
of  a  diagram  of  the  two  patterns  together  which  appears  in  the  Beam  Antenna 
Handbook  (32) •  The  unidirectional  cosine- squared  power  pattern  can  be  expected 
to  provide  a  fair  estimate  of  the  difference  that  can  be  e3q)ected  between 
measurements  made  of  an  event  with  similar  rlometcrs,  one  being  equipped  with  a 
very  narrow  beam  antenna  directed  vertically  upward,  and  the  other  with  the 
three-element  Yagl. 

0 

Letting  G  =*  cos  fl,  ^  is  given  by 

Jr 

a  a^ 

Ap  «  10  log^Q  e“®  1^1 - (1  -  a)j - 

where 


the  exponential  integral. 

In  the  following  table  is  presented  values  of  the  total  vertical  attenuation 
A,  the  attenuation  that  would  be  recorded  by  a  rlometer  equipped  with  an 
(hypothetical)  antenna  having  a  cosine- squared  power  pattern,  A^,  and  the  ratio 
A/Ap  for  selected  values  of  a,  the  Integral  of  the  specific  absorption. 


a 

A(db} 

Ap  (db) 

A/Ap 

0.01 

0.0434 

0.065 

1.50 

0.05 

0.217 

0.323 

1.49 

0.10 

0.434 

0.639 

1.47 

0.50 

2.17 

3.05 

1.4o 

0.75 

3.26 

4.48 

1.38 

1.0 

4.34 

5.88 

1.35 

1.5 

6.51 

8.60 

1.32 

2.0 

8.69 

11.24 

1.29 

3.0 

13.03 

16.36 

1.26 

5.0 

21.71 

26.26 

1.21 

For  events  producing  attenuations  less  than  1  db,  one  could  relate  the  vertical 


to  the  measured  attenuation  by  A  .  *•  1.;/:.  For 

^  ioeasured 


ISA  S  10  db, 

meas.  ' 

10  S  A  <15  db, 

meas .  ' 


meas. 


1.35A, 


A  ~  I.3QA. 
meas.  ^ 


It  may  be  stated  confidently  that  any  large  discrepancies  vhlch  might  exist 
between  the  measured  abcorptlon  and  the  calculated  vertical  absorption  are  not 
due  to  the  receiving  pattern  of  the  rlometer  antenna. 

CHAPTER  K 


Comparison  with  Experimental  Data 

In  order  that  an  euitual  solar  cosmic  ray  event  be  examined  In  light  of  the 
theoretical  resxilts  shown  In  Figs.  23  through  29  and  represented  by  14.  (7«6) 
and  Fig.  30,  there  must  be  available  adequate  latitude  coverage  of  the  cosmic 
radio  noise  absorption  associated  with  the  phenomexia  and  a  measured  cosmic  ray 
spectrum.  Reid  and  Lelnbach  (25)  have  presented  absorption  data  from  rlometers 
situated  at  College,  Barrow  and  IClng  Salmon,  Alaska,  and  Thule,  Greenland,  for 
the  solar  outbursts  occurring  on  iMay  10  and  July  10,  ik,  and  16,  19^9 •  Various 
radio  propagation,  geomagnetic,  and  other  data  related  to  these  same  events, 
but  collected  at  Klruna,  Sweden,  are  given  by  Ortner,  Egeland,  and  Kilqvlst  (29) 
In  their  paper  cited  eeu'ller.  As  far  as  cosmic  radio  noise  data  Is  concerned, 
there  Is  excellent  polar  cap  coverage  for  the  four  great  flares  of  the  year  1959 

On  the  other  band,  a  cosmic  ray  spectrum  derived  from  data  collected  by 
balloon-bonie  counters  has  been  published  for  but  one  of  the  1959  events,  that 
of  July  l4,  vhlch  was  the  second  of  the  July  series  of  three. 

As  reported  by  Winkler,  Bunrsar,  and  Peterson  (l5),  the  sudden  cosnencement 
magnetic  storm  associated  with  the  July  1^  solar  flare  began  at  0602  UT,  July  15 
At  approximately  O9OO  l/T  balloon-borne  counters  of  their  Flight  IQC-lh,  aloft 
at  celling  altitude  over  Mlxmeapolls,  registered  a  very  spectacular  and  rapid 
rise  In  counting  rate.  This  high  intensity  persisted  for  about  five  hours  and 
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then  began  to  decay.  However,  PUglit  IHC-G  wee  launched  at  0004  UT,  July  15, 
and  ascended  during  the  five-hour  period  of  relatively  uniform  high  Intensity. 

Ihe  solar  cosmic  ray  Influx,  therefore,  was  being  monitored  by  one  flight 
(IQC-i4)  at  celling  altitude  while  the  particle  range  spectrum  was  being  obtained 
by  a  subsequent  flight  (IGC-O).  The  cosmic  ray  energy  spectrum  derived  from  the 
data  obtained  through  these  flights  Is  given  by  VJlnklcr,  Bhavsar,  and  Peterson  as 

J(Ej^)  1.05  X  10®E^”^*^  cm"^-sec”^-ster"^  ,  (7*l) 

where  88  <  E^  <  3OO  Mev. 

Absorption  data  taken  at  various  locales  In  the  north  polar  regions  during 
the  July  l4  event  are  shown  in  Fig.  31.  The  data  from  Ft.  Barrow,  niule,  C!ollege, 
and  King  Salmon  were  copied  from  Fig.  6  of  the  paper  of  Reid  and  Leinbach  (2^) 
on  polar  cap  absorption  events.  Ihese  authors  presented  the  absorption  data  as 
superimposed  curves.  &  Fig.  31  of  this  thesis  the  absorption  curves  are  pre¬ 
sented  separately  and  on  a  somewhat  different  scale  than  that  used  by  Beld  and 
leinbach.  The  Klruna  data  were  abstracted  from  Fig.  6  of  the  paper  on  very  low 
frequency  propagation  by  Ortner,  Egeland,  and  Hulqvlst  (29).  lliese  data,  as 
pre  ?nted  in  their  paper,  were  plotted  with  respect  to  Mld-Buropean  Time,  which 
leads  UhlverseQ.  Time  by  one  hour.  Ihls  correction  was  not  attempted  in  Fig.  31, 
the  absorption  data  were  copied  as  presented  in  all  cases  as  faithfully  as 
possible. 

It  Is  evident  that  the  upper  atmosphere  was  still  being  affected  by  a  par¬ 
ticle  bombardment  due  to  the  solar  dlstiixbance  of  July  10  at  O329  UP,  July  l4, 
the  onset  time  of  the  second  solar  event  of  the  July  series.  This  is  indicated 
by  the  presence  of  a  residual  absorption  of  cosmic  radio  noise  of  about  ^  db  at 
all  rlometer  stations  listed  except  King  Saliaon.  Shortly  following  this  second 
great  flare,  a  recrudescence  of  cosmic  radio  noise  absorption  occvirred,  and  at 
about  1^00  UT,  July  l4,  the  cosmic  ray  bosibardment  due  to  this  second  flare  seems 
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to  have  reached  Its  full  development.  This  last  statement  is  based  on  the  fact 
tliat  all  stations  shov  an  increasing  attenuation  of  cosmic  radio  noise  up  to 
about  1^00  UT.  The  College  rioocter  did  record  a  further  sharp  reduction  In 
received  cosmic  radio  noise  power  at  about  2000  UT.  However,  this  appears  to 
have  been  due  to  strong  localised  auroral  activity,  as  It  Is  observed  at  no 
other  station.  The  College  rlometer  record  Indicates  that  an  absorption  In 
excess  of  20  db  occurred  In  the  time  Interval  2300  l/T,  July  l4,  to  about  1200  UT, 
July  13.  It  must  be  pointed  out,  however,  that  such  large  absorption  values  lie 
at  the  very  limit  of  the  dynamic  range  of  the  Instrument.  Hiometer  estimates  of 
cosmic  noise  attenuation  In  excess  of  15  db  are  inaccurate  and  the  actual 
numerical  values  given  should  not  be  taken  literally. 

The  normal  geomagnetic  cutoff  at  the  Thule,  Gbreenland,  Is  the  lowest  of  all 
stations  listed.  Consequently,  the  Thule  rlometer  should  be  expected  to  show  the 
largest  absorption  during  a  solar  proton  event.  Yet  the  composite  record  of  the 
polar  cap  absorption  (Fig.  31)  indicates  that  the  absorption  recorded  at  Thule 
wets  consistently  less  than  that  recorded  at  all  stations  other  than  the  most 
southern  geoaiagnetlcally.  King  Salmon.  For  example,  at  15OO  UF,  July  l4,  the 
following  values  for  cosmic  radio  noise  absorption  were  recorded  at  the  several 
stations: 

Thule,  Greenland  9*5  db 

Ft.  Beurow,  Alaska  16  db 

Klruna,  Sweden  1^  db 

King  Salmon,  Alaska  0  db 

Even  If  one  were  to  attribute  3  to  4  db  absorption  to  auroral  activity  at  the 
intermediate  stations  (Pt.  Barrow,  Klruna,  and  College),  the  whole  Thule  record. 
If  accepted  at  face  value,  remains  an  anomaly.  That  the  accuracy  of  this  record 
should  In  fact  be  held  suspect  is  indicated  by  the  excess  of  absorption  regis¬ 
tered  at  King  Salmon  during  the  magnetic  stozm  (~  12  db)  over  that  recorded  at 
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the  same  time  at  Ihule  (~  10  db).  It  should  be  meutioned  In  this  respect  that 
the  Th\ile  rlotaeter  was  not  under  scientific  surveillance  during  these  July  events^ 
situated  as  It  vas  at  a  rather  remote  outpost,  but  was  checked  from  time  to  time 
by  U.  S.  Azmy  personnel.  During  most  of  the  period  of  Its  operation.  It  vas 
unattended.  It  would  not  be  unreasonable  to  suppose  that  the  Ihule  rlometer  was 
operating  under  a  bias  error  of  at  least  2  db  during  the  July  series.  Okie  may 
compare  the  absorption  record  of  the  Bolar  Blackout  following  the  solar  flare 
of  May  10,  1959»  from  the  Thule  rlometer  with  its  record  of  the  July  events.  Oie 
absorption  at  Thule  during  the  May  event  was,  according  to  the  record,  consis¬ 
tently  greater  than  that  recorded  at  the  other  polar  stations.  It  would  seem 
that  a  systematic  malfunction  of  the  IDiule  rlometer  during  July  Is  a  likely 
eiqplanatlon  of  this  otherwise  unexpUcable  behavior. 

Were  one  to  attempt  to  ascertain  the  cosmic  radio  noise  absorption  due  to 
solar  protons  at  l^OO  UP,  July  1^,  taking  Into  account  auroral  contamination  of 
the  absorption  data  and  a  likely  bias  error  In  the  Thule  rlometer  record,  one 
would  arrive  at  the  following  values  for  solar  cosmic  ray  Induced  cosmic  radio 
noise  absorption  for  the  various  polar  cap  rloamter  stations  listed. 


Station 

Rlometer  Beading 

Correction 

Absorption  due 
to  solar  cosmic 

rays 

Thule 

9.5  db 

2  db  (bias  error) 

11.9  db 

Pt.  Barrow 

16  db 

-3,-4  db  (aurora) 

13,  12  db 

Klruna 

15.5  db 

-3,-4  db  (aurora) 

I2.9>  11.9  db 

College 

l4  db 

-3,-4  db  (aurora) 

11,  10  db 

King  Salmon 

0  db 

— 

0  db 

Is  there  any  corroborative  evidence  that  these  adjusted  absorption  figures 
are  In  fact  correct?  To  find  such  evidence,  one  may  search  the  comoposlte  absozp- 

t 

tlon  record.  Fig.  31,  for  a  time  period  of  minimum  auroral  activity,  but  during 
which  one  could  assume,  with  some  Justification,  that  the  same  solar  proton  flux 
vas  responsible  for  the  cosmic  noise  absorption  recorded  then  as  In  the  adjusted 
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record  In  qviestlOD.  period  shortly  follovrlng  the  tljne  that  the  nwxlmnm 
cosmic  noise  absorption  vas  recorded  at  King  Salioon  seems  to  be  one  of  reduced 
auroral  activity.  All  three  auroral  zone  stations  (College^  Klruna,  and  Pt. 
Barrov)  show  rather  sudden  recoveries  at  about  13OO  July  15.  abruptness 
of  these  recoveries  and  the  leu^  of  change  In  the  corresponding  niule  record 
Indicate  a  diminution  of  auroral  activity  rather  than  a  change  in  the  solar 
cosmic  ray  flux.  The  cheusges  that  had  occurred  at  College^  Klruna  and  Pt. 

Barrow  at  this  time  are  to  be  contrasted  with  those  evidenced  at  all  stations 
somewhat  later  (commencing  at  about  16OO  IKT,  July  15)  which  were  much  more 
systematic  end  regular,  l^ese  latter  changes  seem  to  represent  a  recovery  due 
to  a  decay  in  the  solar  cosmic  ray  Intensity  rather  than  a  decrease  In  aurorol 
activity. 

At  1300  UT,  July  13 i  about  two  hours  following  the  sudden  recoveries 

recorded  at  College,  Klruna  and  Pt.  Barrow  Just  alluded  to,  the  followlxig  values 

of  cosmic  radio  noise  absorption  were  recorded  at  the  various  riometer  stations: 

Thule  10,5  +  2  (bias  error)  »  12.5  * 

Pt.  Barrow  10  db 

Klruna  11  db 

College  10  db 

King  Salmon  9  db 

The  cosmic  radio  noise  absozptlon  values  presented  above  show  consistency 
with  the  adjusted  values  of  I5OO  IfF,  July  1^,  appearing  earlier  In  this  section, 
but  with  King  Salmon  showing  9  db  cosmic  noise  absorption.  This  latter  reading 
Is  Indicative  of  a  breakdown  of  the  geomagnetic  cutolf  that  existed  before  the 
onset  of  the  magnetic  stozm  associated  with  the  July  l4  flare. 

In  sxmmory,  the  above  observations  seem  to  support  the  following  assumptions 
concerning  the  solar  event  of  July  1^,  1959: 

(1)  The  fully  developed  solar  proton  Influx  bombarding  the  upper  atmosphere 


during  the  time  Interval  I5OO  UP,  July  lU,  to  about  18OO  UP,  July  15,  1959>  can 
he  represented  hy  the  Integral  spectrum  obtained  by  VJlnkler,  Shavsar,  emd 
Beterson  from  balloon-borne  counters  at  Minneapolis,  Minnesota,  shortly  after 
the  sudden  comnencement  of  geomagnetic  storm  of  July  15,  19^9 •  ^la  spectrum 
can  be  presented  In  the  fom  of  a  power  law  as  follows: 

J(E)  =  1.05  X  10^^"^*^  protons/cm^-sec-ster,  (7-l) 

where  68  <  <  3OO  Mev. 

(2)  Accotintlng  as  far  as  possible  for  absorption  Induced  by  auroral  activ¬ 
ity  during  the  solar  cosmic  ray  event  and  Inaccuracies  In  the  rlometer  data, 

one  woiild  conclude  that  the  cosmic  radio  noise  absorption  due  to  the  solar  proton 
bombardment  during  the  height  of  the  July  ik  event  was  virtually  constant  over 
the  polar  cap,  and  heui  a  numerical  magnitude  of  10  to  12  db. 

(3)  ^  geomagnetic  cutoff  was  operative  from  the  onset  time  of  the  July  l4 
event  to  the  time  of  the  sudden  commencement  of  the  geomagnetic  storm  of  July  1^. 
This  Is  directly  adduclble  from  the  cosmic  radio  noise  absorption  record  from 
King  Salmon,  Alaska,  which  shows  no  cosmic  radio  noise  absorption  up  to  the  time 
of  the  S.C. 

It  Is  Immediately  obvious  that  the  theoretical  variation  of  cosmic  radio 
noise  absorption  with  geomagnetic  latitude  (Fig.  29)  Is  not  reflected  In  the 
absorption  data  of  this  event  even  when  allowances  are  made  for  auroral  con¬ 
tamination  and  rlometer  malfunction.  It  must  be  remembered  In  this  respect  that 
the  theoretical  variation  of  cosmic  radio  noise  absorption  as  determined  In  this 
thesis  Is  predicated  on  the  assumption  that  extrapolation  of  the  solar  cosmic 
ray  spectrum  to  low-energy  values  is  permissible.  The  lack  of  variation  of 
cosmic  radio  noise  absorption  over  the  polar  cap  in  the  presence  of  a  geomag¬ 
netic  cutoff  Is  direct  evidence  of  the  Inadmissibility  of  this  assunqptlon  In 
this  case  at  least.  The  solar  proton  Integral  spectrum  must  have  begun  to 


flatten  somevhere  In  the  nelSiihorhood  of  90  or  100  Mev  to  produce  the  distribu¬ 
tion  of  polar  cap  absorption  Inferred  from  the  data.  Ihls  statement  may  seem  to 
conflict  vlth  the  claim  of  Winkler  et  al.  that  their  power- lav  spectrum  Is  valid 
to  88  Mev.  However,  It  Is  the  only  explanation  that  accounts  for  the  observed 
variation  of  cosmic  noise  absorption  In  the  presence  of  the  nozmal  (Quenby- 
Webber)  cosmic  ray  cutoff.  Were  one  to  suppose  that  the  geomagnetic  cutoff 
actually  operative  on  July  l4,  1959>  vas  quite  different  from  that  normally 
Imposed  by  a  dipole- like  geomagnetic  field,  he  would  find  that  this  would  not 
obviate  the  theoretical  necessity  for  a  flattening  of  the  cosmic  ray  spectrum. 

It  would,  however,  permit  a  lower  estimate  to  be  made  of  the  energy  value  at 
\dilch  the  flattening  of  the  spectrum  began.  Unfortunately  there  is  no  direct 
evidence  that  the  cosmic  ray  cutoff  was  anything  other  than  that  normally  opera¬ 
tive  over  the  polar  cap.  For  protons  the  cutoff  at  College  Is  112  Mev.  That 
there  had  been  virtually  the  same  cosmic  noise  absorption  observed  at  Ft.  Barrow 
and  Thule  as  that  recorded  at  College  requires  a  flattening  of  the  solar  proton 
Integral  spectrum  at  energy  values  at  least  comparable  to  the  geomagnetic  cutoff 
for  protons  at  College. 

Can  a  correct  estimate  of  the  cosmic  noise  absorptlca  during  the  event  at 
a  given  station  be  made  from  the  results  of  calculation,  based  as  they  are  on 
the  values  of  the  various  recoablnatlon  rate  coefficients  adopted  In  this  thesis? 
To  answer  this  question,  the  solar  proton  Integral  spectrum  obtained  by  Winkler, 
Bhavsar  and  Beterson  can  be  \ised  with  the  cutoff  for  protons  obtained  from  the 
theory  of  Quenby  and  Webber  In  Eq.  (j.S)  to  obtain  the  total  vertical  attenuation 
at  College,  Alaska,  at  1500  UT,  July  l4.  The  conposlte  absorption  record.  Fig. 
31,  Indicates  that  the  solar  proton  bombardment  vas  fully  developed  at  that  time 
and  the  King  Salmon  record  shows  the  existence  of  a  geomagnetic  low-energy  cut¬ 
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Substituting  the  proper  numerical  values  for  J(ll2  Mev),  7,  and  E  into 

c 

(7>6);  one  obtains: 

J(112  Mev)  -  1.05  X  =  119.7  , 

r  =  1  +  2.9  =  3.9  , 


and 


Eg  =  112  Mev  , 

M(Eg,7)  =  0.80  (from  interpolation  in  Pig.  3O)  > 

log^Q(A^(ll2)3*^/l35.7]  =  2(3.9  -  1)  +  0.80  , 


where  A  is  the  total  vertical  absorption  in  db.  This  yields 


A  =  2.18  db  . 

The  cosmic  radio  noise  absorption  recorded  by  the  College  riometer  at  I5OO 
July  l4>  was  l4  db.  If  one  takes  into  account  the  absorption  due  to  auroral 
activity,  as  was  suggested  earlier  in  this  section,  he  must  subtract  3  to  4  db 
from  the  vedue  l4  db  given  in  the  record.  A  further  correction  for  the  finite 
beam  width  of  the  antenna,  discussed  in  an  earlier  section,  reduces  the  absorp¬ 
tion  figure  by  about  2  more  db.  Obis  leaves  a  total  vertical  absorption  of 
about  9  db,  or  a  discrepancy  of  about  7  db  between  the  actual  total  vertical 
attenuation  of  cosmic  radio  noise  that  can  be  directly  attributed  to  the  Influx 
of  solar  protons  during  the  July  l4  event  and  that  deduced  from  the  cosmic  ray 
spectrum  measured  in  the  course  of  the  event.  That  less  attenuation  is  pre¬ 
dicted  by  a  calculation  from  the  cosmic  xay  spectrum  than  is  actually  obsezved 
is  due  either  to  the  use  of  an  Incorrect  electron-neutral  collision  frequency 
profile  in  the  calculation,  or  to  an  underestimation  of  the  equilibrium  electron 
density.  Tbe  second  alternative  is  almost  certainly  the  case,  since  the  numeri¬ 
cal  determination  of  the  electron  density  due  to  a  solar  proton  bosibardfflBnt 
depends  rather  critically  on  the  values  adopted  for  the  various  atmosphere  rate 
coefficients.  It  was  indicated  in  an  earlier  section  of  this  thesis  that 
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\incert«lntles  of  an  order  of  magnitude  or  more  can  exlat  In  the  valuea  one  ml^ht 
choose  to  assign  to  these  coefficients. 

The  curves  representing  the  specific  attenuation  given  in  Fig.  10  show  that 
the  atmospheric  height  range  of  greatest  importance  for  the  absorption  of  cosmic 
radio  noise  is  that  of  30  to  63  km.  At  these  heij^ts  the  most  important  reaction 
occurring  to  establish  the  ultimate  value  of  the  electron  density  caused  by  the 
production  of  a  free  electron  by  solar  cosmic  rays  la  the  attachment  of  electrons 
to  oxygen  molecules.  The  rate  coefficient  of  this  three>body  process  has 
recently  been  revised  downward  from  a  value  of  approximately  2  x  lO'^nCOg)— 
[nCOg)  is  the  nvudber  density  of  molecular  oxygen]  used  by  Bailey^  to  that  of 
1.3  X  10*^  n(0g),  quoted  by  Nicolmt  and  adopted  for  use  in  the  present  woxic.  The 
results  of  the  investigation  at  cosmic  radio  noise  attenuation  caused  by  solar 
cosmic  rays  presmsted  In  this  thesis^  indicate  that  a  xeinvestlgatlon  of  this 
electron  attachment  process  is  needed.  A  further  downward  revision  of  the  value 
of  the  rate  coefficient  of  this  process  would  bring  agreement  between  the  cal¬ 
culated  and  observed  cosmic  radio  noise  absorption  due  to  a  solar  corpuscular 
bOid>ardnent  of  the  earth's  atmosphere. 

CHAFIBR  X 

Suamry  and  Conclusioos 

In  the  course  of  this  thesis  there  has  been  described  a  method  of  determin¬ 
ing  the  changes  produced  in  the  upper  atmosphere  by  an  influx  of  solar  cosmic 
rays.  Such  a  boeibardment  brings  into  being  a  condition  of  abnormal  ionisation 
in  a  region  below  the  normally  occurring  D  layer.  This  is  directly  dlsoemable 
as  a  change  in  the  radlo-vave  propegatlon  characteristics  of  the  upper  atmos- 
phere^  the  region  affected  beccadng  one  of  relatively  strong  nondeviatlve  absorp* 
tion.  This  absorption  is  confined  to  the  earth's  polar  regions.  A  quantitative 
measure  of  the  absorbing  power  of  the  upper  atmosphere  under  the  Influence  of  a 


solar  cosmic  ray  bombardment  can  be  obtained  throng  continuous  monitoring  of 
the  Intensity  of  the  galactic  X'adlo  noise  reaching  the  earth's  surface  by 
special  receiving  systems  (rlometers)  distributed  over  the  geographical  area 
affected. 

The  Intensity  and  distribution  of  cosmic  radio  noise  absorption  Is  deter¬ 
mined  theoretically  In  this  thesis  under  the  following  basic  assumptions: 

(1)  The  solar  cosmic  ray  differential  energy  spectrum  can  be  represented 
by  a  power  lav  of  the  form 

N  -  , 

where  a  and  7  are  constants  and  is  the  free-space  particle  energy.  (This 
Implies  an  integral  spectrum  of  similar  form  which  is  consistent  with  results 
derived  from  balloon  measxirements ) . 

(2)  A  low-energy  cutoff,  E  ,  is  Imposed  on  the  spectrum  by  the  earth's 

c 

magnetic  field.  The  value  of  the  cutoff  energy  for  a  given  pair  of  geomagnetic 
coordinates  Is  obtained  from  the  theory  of  Qnenby  and  V7ebber. 

(3)  theoretical  electron  density  profile  existing  In  the  region  of 
abnonsal  Ionization  can  be  deduced  from  the  assui^ption  that  a  state  of  quasl- 
equlllbrlum  exists  between  the  production  and  recombination  of  free  electrons. 
The  Individual  electron  production  and  recosibinatlon  processes  assumed  to  be  of 
IsqportaQce  In  the  establishment  of  the  electron  density  profile  are  dissociative 
recariblnatlon  and  electron  attachment  to  molecular  oxygen  (loss  mechanisms) 

and  ptaotodetachment  and  coUlslonal  detachment  of  electrons  from  heavy  negative 
Ions  (secondary  production  mechanisms).  A  process  of  mutual  neutralization 
between  positive  and  negative  Ions  also  occxxrs.  This  process  Is  important  In 
establishing  the  negative- lon-to-electron  ratio. 

Calculations  of  the  total  vertical  absorption  of  cosmic  radio  noise  were 
pexfonaed  using  three  trial  solar  proton  spectra  for  four  low-energy  cutoff 
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values  used  correspond  to  the  geomagnetic  cutoffs  at  Pt.  Barrow,  Ft.  Yukon, 
College,  and  Farewell,  Alaska,  respectively.  The  results  of  the  calculations 
were  sufficiently  well  ordered  to  permit  the  construction  of  ahsorptlon  curves 
which  are  given  in  Figs.  23  through  28.  flirthermore,  it  was  possible  to  obtain 
an  en^irlcal  formula  of  simple  form  relating  the  total  vertical  absorption  of 
cosmic  radio  noise  to  the  solar  cosmic  ray  spectrum  and  the  low-energy  cutoff. 
Ibis  formula  represents  a  continuation  of  the  toted,  vertical  absorption  for 
values  of  y  and  Ec  other  than  those  used  orlglxrally. 

The  results  obtained  from  the  theory  are  compared  with  observations  made  on 
the  solar  cosmic  ray  event  of  July  l4,  19^9 •  After  reasonable  8d.lowances  are 
made  for  auroral  effects  and  possible  rlometer  malfunction  during  the  event,  the 
following  discrepancies  between  the  theoretical  results  and  observation  are  shovm 
to  exist: 

(1)  The  distribution  of  cosmic  radio  noise  absorption  over  the  polar  cep 
did  not  conform  in  this  event  with  that  derived  under  the  assumptions  given 
above.  The  absorption  observed  was  virtually  constant  while  the  cadculated 
distribution  of  absorption  indicates  that  a  strong  variation  with  geomagnetic 
latitude  should  exist. 

(2)  The  theory  yields  considerably  less  absorption  for  a  known  solar 
cosmic  ray  flux  at  a  given  location  than  that  actuedly  observed. 

The  discrepancy  in  latitude  variation  is  indicative  of  a  failure  of  the 
power- law  speotrum  to  hold  for  solar  proton  energies  much  below  100  Mev  in  this 
particular  event.  This  means  that  extrapolations  of  the  spectrum  to  energy 
values  less  than  those  directly  measurable  can  be  an  invalid  procedure.  A  lack 
of  variation  of  cosmic  radio  noise  absorption  with  geomagnetic  latitude  implies 
a  flattening  of  the  solar  cosmic  ray  spectrum  in  the  neighborhood  of  the  hipest 
cutoff  energy  displaying  the  constant  absorption. 

The  failure  of  the  theory  to  account  for  the  magnitude  of  the  cosmic  radio 
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noise  absorption  obser/ed  must  be  attributed  to  one  or  more  Incorrect  values  of 
the  atmospheric  recombination  rate  coefficients.  The  particular  coefficient 
most  suspect  Is  that  of  the  three-body  electron  attaclunent  process  to  molecular 
o^^en. 

nils  process  la  the  most  Important  one  occurring  in  the  atmospheric  region 
of  greatest  absorption;  therefore,  an  erroneous  value  of  Its  rate  coefficient 
vould  Introduce  an  appreciable  error  in  the  final  results  of  the  absorption 
calculation.  Furthermore,  this  rate  coefficient  has  recently  undergone  an 
order-of -magnitude  revision  in  Its  numerical  value.  A  large  degree  of  uncer¬ 
tainty  still  exists  In  the  numerical  value  one  should  ascribe  to  this  quantity. 

APPENDIX 

z  h 


The  particle  whose  trajectories  and  atmospheric  stopping  helots  were 
detenslned  In  the  analysis  to  be  described  were  considered  to  be  protons  having 
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initial  energies,  in  tlie  range 

20  <  E^  <  260  Mev  . 

Belativlstic  effects  vere  neglected  and  the  classical  equations  of  motion  were 
used  throughout. 

p,  q,  h  be  a  cartesian  coordinate  system  sltmted  on  the  earth  at  a 
given  geomagnetic  latitude;  h  la  positive  upvard  In  the  direction  of  the  loccQ. 
zenith.  This  system  will  be  called  the  earth  coordinate  syaten. 

I«t  X,  y,  z  be  a  cartesian  coordinate  system  so  situated  that  z  Is  parallel 
but  oppositely  directed  to  the  local  terrestrial  magnetic  field  vhieh  Is  con¬ 
sidered  (in  the  small)  to  be  of  constamt  magnitude  and  direction,  ^is  system 
vlU  be  called  the  wwjmetle  field  coordinate  system.  Ihe  axes  Oh,  Oz,  and  Op 
are  coplaner,  and  X  Is  the  angle  between  the  zenith  direction  (h)  and  the  direc¬ 
tion  of  the  z  axis  (opposite  to  the  earth's  field  direction).  !Qie  q  and  y 
coordinates  for  the  two  systems,  respectively,  are  coUinear. 

Ihe  cylindrical  coordinate  system.  In  whose  terms  the  equations  of  motion 
of  a  proton  of  mass  m  and  charge  e  traveling  through  the  atmosphere  are 
e3qi>ressed.  Is  related  to  the  magnetic  field  coordinate  system  throu^  the 
following  equations. 

X  -  R  cos  9  , 

y  »  R  sin  6  ,  (A.l) 

z  «  z  . 

Ihe  equations  of  motion  themselves  are 

mR  -  mRR^  ♦  eflRB^  “  ' 

aB9  +  2bB9  -  eRB^  “  ^0  *  (A. 2) 

m  *>  Q  , 
z 

where  m  Is  the  mass  and  e  the  charge  of  the  particle,  B^  Is  the  earth's  field 
and  If  Is  a  dissipative  force  given  In  terms  of  the  energy- loss  function  for 
protons  (Fig.  1,  Curve  A).  !Ihe  dot  denotes  differentiation  with  respect  to  the 
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time,  nie  components  of  ^ ere  given  as  follows: 


d£ 

Q  =  -  p  —  ? 
®  dS  ® 


dE 

A 

-  P  —  e 


dS 

dE 


e 


Q  B  -  p  —  e 
*  dS  * 


®R  “  “  “"TT  > 
®  dS  V 


.  <113  V, 

e*,  -  -  P  —  7-7  > 
®  dG  V 


dE  V 

6(5  =  -P  —  77-  > 
dS  V 


(A.3) 


vhcre  P  is  the  air  density  at  a  height  hj  dE/dS  is  the  energy- loss  function 
(with  a  small  change  in  notation)  given  in  Fig.  1,  Curve  A;  e^,  e^,  e^  are  unit 
vectors  pointing  in  the  positive  radial> azimuthal  and  axial  directions  in  the 
magnetically  centered  cylindrical  coordinate  system;  eg  ia  a  imit  vector  point¬ 
ing  in  the  instantaneous  direction  of  the  particle  motion;  v^  v^,  v^  are  the 
cosqponents  of  the  particle  velocity  in  the  R>  9,  and  z  diz«ctlons,  and  |v|  is 
the  absolute  value  of  the  total  velocity. 

The  following  substitutions  are  made  for  the  dependent  variables  R  and  0 
and  for  the  Independent  variable^  t  (the  time), 

B  -  .  t  , 

9  ■  9  +  OD^t  ,  (a. 4) 

T  »  O^t  . 

R^  is  the  initial  radius  of  the  circular  spiral  along  vhich  the  proton  is  moving 

before  being  slowed  down  in  the  atmosphere,  and  oi^  is  the  initial  gyrofrequency. 

OD^  »  eB^m  .  (A. 5) 

Making  these  substitutions  into  the  equations  of  motion  (Bq.(A.2)],  one 

obtains  the  following,  after  a  certain  amount  of  reduction  and  siatpllflcatlon: 

.2« 


dt  /dq)\/dq)  \  d( 

— P  -  (R  +  U  —  )(  —  +  1)  -  -  P  —  , 

dT^  °  A  dT  /  N  dT  /  dT 
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(2(ci9/dT)  +  1]  di/ix 


9 


(A.6) 


vhere 

Qp(h)(dJ;/dS)(E) 

p  ■  ■  ^  (A«7) 

iiBfl^((dt/dT)2  +  Rq  +  O^Udtp/dT)  +  1]^  +  (dz/dT)^)^^  * 

In  the  actual  numerical  vork,  the  following  units  were  used:  length  in 

kilometers,  p(h},  the  density,  in  graas/cm^,  m,  the  mass,  in  grams,  e,  the  charge, 

in  couloofbs,  0),  the  gyrofre4,uency,  in  radlans/sec,  (dE/dS}(E),  the  energy  loss, 

in  Mev  per  ga/cn?,  a  is  a  numerical  constant  whose  value  depends  on  the  units 

used,  (it  was  assigned  a  value  of  1.60  x  10~^  which  is  consistent  with  the 

units  employed  in  the  calculation.)  t,  dimensionless,  F,  dimensionless. 

'Dene  initial  conditions  Imposed  on  the  equations  of  motion  (Eq.  A.6)  to 

deteznlne  a  trajectory  are:  h^  >  the  initial  (starting  height  of  the  proton 

ahove  the  earth;  p^,  q^,  h^  >  the  components  of  the  particle  velocity  at  height 

h^  in  the  earth  coordinate  system.  These  quantities  along  with  a  value  for  the 

angle  X  are  sufficient  to  specify  the  trajectory.  The  actxial  physical  location 

of  the  trajectory  is  irrelevant.  In  the  magnetic  cylindrical  coordinate  system 

the  initial  conditions  on  R.,  R.,  0.,  and  9.  are:  R.  ■  R^  (R.  *  the  ixiitial 

ill  i  101 

spiral  radius);  R^  «  0  (initially  the  particle  trajectory  is  a  clreiilar  spiral); 

9^  •  (t  »  0);  (the  free-space  gyrofrequency). 

In  terms  of  the  transfonaed  varlahles  we  have 


dl 

dx 


“  0  , 

-  0  , 

1 


dp 

— •  ■  0  . 

dx  1 


(A.8) 
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Bquatlono  (A. 8)  hold  for  any  trajectory.  A  particular  trajectory  la  deter¬ 
mined  by  the  Initial  values  of  z,  z,  and  (p,  idilch  are  derived  from  the 
initial  values  h^,  p^,  4^^  h^. 

The  initial  value  of  z  can  be  expressed  in  terms  of  and  <p^. 

z^  ■  R^  cos  q>^  tan  ^  +  h^^  sec  X  .  (A. 9) 

NOW 

®o  -  > 

where  ^ol  is  the  kinetic  energy  associated  with  the  motion  tremsverse  to  the 

magnetic  field  at  heigdit  h^,  and  is  the  component  of  the  initial  velocity 

(the  velocity  at  height  h^)  parallel  to  the  magnetic  field.  After  a  certain 

amount  of  algebraic  manipulation,  it  can  be  shown  that 

R^  =  (l/c^)[(Pj^  cos  X  -  hj^  sin  X)^  +  .  (A. 10) 

Ibtis  Bq,.  (a.  10)  expressed  R^  in  terms  of  initial  values  of  the  veloci'ty  in  the 

earth  coordinate  system. 


so 


dz 

dt 


z  “  p  sin  X  +  h  cos  X  , 


=  —  (p.  sin  X  +  h.  cos  X)  , 
T»0  <0^ 


(A.U) 


and  finally,  because  the  initial  velocity  vector  must  be  normal  to  the  initial 


radius  vector,  R^,  ve  have 


-  tan‘^(-Xj^/yj^)  , 

<Pj^  «  tan"^l(hj^  sin  X  -  p^^  cos  X)/4^]  .  (A.12) 

Care  must  be  exercised  in  chooslxig  the  correct  brcmch  of  this  function  [Eq. 

(a. 12)1.  Thus  all  initial  quantities  of  importance  are  given  in  terns  of  the 
necessary  initial  conditions  in  the  earth  coordinate  system. 

To  detemlne  a  trajectory,  one  must  solve  the  equations  of  motion,  Eqs. 

(a. 6),  numericadly  with  the  above  initial  conditions.  The  changes  of  variable 
in  the  differential  equations,  especially  in  the  independent  variable,  were  made 


65 


primarily  for  scaling  purposes  to  facilitate  the  nianerlcal  integration*  Final 
values  computed  were  and  for  stopping  particles,  l.e.,  for  the  con¬ 

dition  that  the  particles'  kinetic  energy  equals  zero.  A  further  transformat Icm 


of  coordinates  gives  the  stoppliag  height  h^, 

h^  »  Zj  cos  X  -  Xj  sin  X  . 

The  contents  of  this  Appendix  can  be  summarized  as  follows: 


(A. 13) 


dh 

— 9  -  (R^  +  i) 

dT  dT  X  dT 


-  F 


d^cp  [2(d9/dT)  +  1]  dg/dT 
— -  + -  *  _  F 


dg 

dT 

dep 


dT 


d^z 


(rO' 


(A.6) 


dT 


dz 

F  —  ; 

dT 


1.60  X  10“^^p(h)(dE/dS)(E) 


,^{(dg/dT)^  +  (R^  +  g)^[(d(p/dT)  +  l)^  + 


(A.7) 


moig 

The  particle  kinetic  energy  in  Mev  is  given  by 
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E(M.v)  .  6.E5  X  10«  :i[(  jy  .  (E„  -  S)=(  *  ij  • 


dep 

dT 


dz  \2 

dT 


Initially 


0  I  — 


dg 

dT 


0  ,  — 
dT 


0  , 


(A.8) 


and  in  terms  of  the  initial  values  h^,  p^,  q^^^  the  other  necessary  initial 
values  are  given  as  follows: 

Rjj  =  (iAu^)l(Pj^  cos  X  -  hj^  sin  X)^  +  , 


=  R^  cos  qj^tan  X  +  hj^  sec  X  , 


dz 

dT 


,  e  a  . 

«  —  (p.  sin  X  +  h .  COB  X)  , 

»  tan”^[(hj^  sin  X  -  p^^  cos  X)/q^]  , 


(A.IO) 

(A.9) 

(A.n) 


(A.12) 
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13ia  coordinate  transformations  'ised  were 

X  ■  p  cos  X  -  h  sin  X  , 

y  =  4  ,  (a. 15^ 

z  =  p  sin  X  +  h  cos  X  , 

and  Its  Inverse.  Also,  the  cylindrical  system  Is  related  to  the  cartesian  mag¬ 
netic  system  through 

X  “  R  cos  9  , 

y  =  R  sin  9  ,  (A.l) 

z  B  z  , 
and 

hj  »  cos  X  -  x^  sin  X  .  (A.13) 

In  addition  to  the  trajectory  and  stopping  height  calculations,  a  particle 
range  calculation  vas  made, 

S  y  ^  S 

Ibis  vas  done  to  provide  a  check  on  the  correctness  of  the  numerical  work.  Ibe 
range  In  air  of  a  particle  vlth  Initial  energy  Is  known  and  the  confpited  range 
can  be  directly  compared  vlth  known  values . 

The  equations  of  taotlon  and  Initial  conditions  used  in  the  forms  given  above 
allowed  one  to  study  the  structure  of  particle  trajectories  as  well  as  to  deter¬ 
mine  particle  stopping. heights.  Of  considerable  Interest  was  the  amount  of 
spiraling  the  particle  performed  in  Its  Journey  down  Uirough  the  atmosphere.  The 
Initial  helc^t  (h^)  chosen  to  begin  the  numerical  Integration  vas  10$  km.  The 
initial  velocities  were  chosen  in  accordance  vlth  various  stipulated  Initial 
directions  vlth  respect  to  the  earth's  field  and  initial  kinetic  energies.  The 
numerical  resiats  displayed  negligible  proton  energy  loss  occurring  at  heights 
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greater  than  about  65  Im,  but  at  heights  less  than  this  the  rate  of  loss 
becomes  great,  the  particles  being  stopped  rather  abruptly.  Ibe  nvanerlcal 
solutions  also  indicate  that  for  solar  cosmic  rays  there  Is  no  spiraling  in  the 
atmosphere— the  rectilinear  approximation  to  the  actual  particle  trajectory  is 
almost  always  adequate. 


REFERENCES 


68. 


1.  P.  Meyer,  E.  Parker,  and  J.  A.  Simpson,  Phys.  Rev.  104,  768 
(1956). 

2.  B.  Peters,  J.  Oeophys.  Res.  155,  (1959). 

5.  E.  P.  Ney,  Ann.  Rev.  Nucl.  Science  10,  46l  (i960). 

4.  I.Q.Y.  Bulletin,  Natl.  Acad.  Sciences,  No.  35,  May,  I960. 

5.  A.  P.  Mitra  and  C.  A.  Shain,  J.  Atmos.  Terrest.  Phys.  4,  204 
(1955). 

6.  M.  Nlcolet  and  A.  C.  Alkin,  J.  Oeophys.  Res.  1469  (i960). 

7.  B.  Rossi,  High  Energy  Particles,  (Prentice-Hall,  Inc.,  New 

York,  1952). 

8.  R.  L.  P.  Boyd  and  M.  J.  Seaton,  eds..  Rocket  Explorations  of 
the  ^pper  Atmosphere  (Pergamon  Press  Ltd.,  London,  1954). 

9.  N.  Nlcolet,  Phys.  Fluids  2,  95  (1959). 

10.  Rocket  Panel,  Phys.  Rev.  1027  (1952). 

11.  E.  Segr^,  ed..  Experimental  Nuclear  Physics  (John  Wiley  and 
Son,  New  York,  1955)  Vol.  I. 

12.  D.  K.  Bailey,  Proc.  I.R.E.  255  (1959). 

15,  R.  R.  Brown  and  R.  0.  D'Arcy,  Phys.  Rev.  Letters  2*  590 
(1959). 

14.  K.  A.  Anderson,  R.  Amoldy,  R.  Hoffman,  L.  Peterson,  and 
J.  R.  Wlnckler,  J.  Oeophys.  Res.  64,  II55  (1959). 

15.  J.  R.  Wlnckler,  P.  D.  Bhavsar,  and  L.  Peterson,  J.  Oeophys. 
Res.  66,  995  (196l). 

16.  J.  J.  Quenby  and  W.  R.  Webber,  Phil.  Mag.  4,  90  (1959). 

17.  J.  E.  Kasper,  Suppl.  Nuovo  Clmento  11,  1  (1959). 


69. 

18.  E.  P.  Ney,  Proc.  International  Conference  on  Cosmic  Ray 
Physics,  Kyoto,  Japan  (196I)  (in  press). 

19.  C.  M.  Crane,  J.  Oeophys.  Res.  66,  111?  (196I). 

20.  A.  V.  Phelps  and  J.  L.  Pack,  Phys.  Rev.  Letters  6,  111 

(1961). 

21.  L.  M.  Channin,  A.  V.  Phelps,  and  M.  A.  Biondi,  Phys.  Rev. 
Letters  2,  344  (1959). 

22.  S.  J.  Smith,  D.  S.  Burch,  and  L.  M.  Branscomb,  Ann.  Oeophys. 

225  (1958). 

23.  H.  Leinbach,  Private  communication. 

24.  D.  R.  Bates  and  H.  S.  W.  Massey,  Proc.  Roy.  Soc.  (London) 
AI87.  261  (1946). 

25.  0.  C.  Reid  and  H.  Leinbach,  Report,  Geophysical  Institute, 
University  of  Alaska,  April,  i960. 

26.  J.  A,  Kane,  J.  Oeophys.  Res.  64,  I33  (1959). 

27.  J.  M.  Anderson  and  L.  Goldstein,  Phys,  Rev,  102.  388  (1956). 

28.  Geomagnetic  Survey  Maps,  Nos.  17OI,  1703,  U.  S.  Hydrographic 
Office  (1955). 

29.  J.  Ortner,  A.  Egeland,  and  B.  Hultqvlst,  I.R.E.  Trans. 
Antennas  and  Propagation  AP-8,  621  (i960). 

30.  B.  Hultqvlst  and  J.  Ortner,  Planetary  and  Space  Science  1, 
193  (1959). 

31.  C.  0.  Uttle  and  H.  Leinbach,  Proc.  I.R.E.  42,  315  (1959). 

32.  W.  I.  Orr,  Beam  Antenna  Handbook,  Radio  Publications  Inc., 


Danbury  Road,  Wilton,  Conn. 


FiaURE  1£0ENDS 


70. 


Pig.  1.  Proton  Energy  Loss  in  Air  vs  Proton  Energy 

Pig.  2.  Electron  Production  Rate  vs  Atmospheric  Height  for  Three 
Solar  Proton  Spectra. 

Fig.  3.  Electron  Production  Rate  vs  Atmospheric  Height  for  Solar 
Protons  and  Alpha  Particles. 

Pig.  4.  Equilibrium  Electron  Density  vs  Atmospheric  Height  under 
Daytime  and  Nighttime  Conditions. 
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Recombination  Coefficient  vs  Atmospheric  Height  under  Daytime  Condi¬ 
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tions. 
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under  Daytime  and  Nighttime  Conditions. 
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Fig.  14.  Constant  Height  Contours  of  Electron  Production  Rate  vs 
Solar  Proton  Cutoff  Energy  for  7-5  Solar  Cosmic  Ray  Spectrum. 

Pig.  15.  Constant  Height  Contours  of  Electron  Production  Rate  vs 
Solar  Proton  Cutoff  Energy  for  7*5  Solar  Cosmic  Ray  Spectrum. 

Fig.  16.  Constant  Height  Contours  of  Electron  Production  Rate  vs 
Solar  Proton  Cutoff  Energy  for  7  ■  7  Solar  Cosmic  Ray  Spectrum. 

Pig.  17.  Constant  Height  Contours  of  Equilibrium  Electron  Densi¬ 
ty  (Daytime)  vs  Solar  Proton  Cutoff  Energy  for  7  ■  3  Solar  Cosmic  Ray 
Spectrum. 

Pig.  18.  Constant  Height  Contours  of  Equilibrium  Electron  Densi¬ 
ty  (Daytime)  vs  Solar  Proton  Cutoff  Energy  for  7-5  Solar  Cosmic  Ray 
Spectrum. 

Pig.  19.  Constant  Height  Contours  of  Equilibrium  Electron  Densi¬ 
ty  (Daytime)  vs  Solar  Proton  Cutoff  Energy  for  7-7  Solar  Cosmic  Ray 
Spectrum. 

Pig.  20.  Constant  Height  Contours  of  Specific  Attenuation  (Day¬ 
time)  vs  Solar  Proton  Cutoff  Energy  for  7-3  Solar  Cosmic  Ray 
Spectrum. 

Pig.  21.  Constant  Height  Contours  of  Specific  Attenuation  (Day¬ 
time)  vs  Solar  Proton  Cutoff  Energy  for  7*3  Solar  Cosmic  Ray 
Spectrum. 

Pig.  22.  Constant  Height  Contours  of  Specific  Attenuation  (Day¬ 
time)  vs  Solar  Proton  Cutoff  Energy  for  7  ■  7  Solar  Cosmic  Ray 
Spectrum 

Pig.  23.  Total  Attenuation  vs  Solar  Proton  Cutoff  Energy  (Day¬ 
time),  7  "  5. 

Pig.  24.  Total  Attenuation  vs  Solar  Proton  Cutoff  Energy  (Day¬ 
time),  7  ■  5. 
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Pig.  25.  Total  Attenuation  vs  Solar  Proton  Cutoff  Energy  (Day¬ 
time),  7  “  7. 

Pig,  26,  Total  Attenuation  vs  Solar  Proton  Cutoff  Energy  (Night¬ 
time),  7  "  3. 

Pig.  27.  Total  Attenuation  vs  Solar  Proton  Cutoff  Energy  (Night¬ 
time),  7*5. 

Pig.  28.  Total  Attenuation  vs  Solar  Proton  Cutoff  Energy  (Night¬ 
time),  7  ■  7. 

Pig.  29.  Total  Attenuation  vs  Geomagnetic  Latitude  (Daytime  and 
Nighttime)  for  Three  Solar  Cosmic  Ray  Spectra. 

Pig.  50.  M(E^,7)  vs  Solar  Proton  Cutoff  Energy  (see  text), 
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Pig.  31.  Cosmic  Radio  Noise  Absorption  Data  from  the  Event  of 
July  14,  1959. 
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